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A MODULAR RADIATIVE TRANSFER PROGRAM FOR GAS 
FILTER CORRELATION RADIOMETRY 

By 

Joseph C. Casas ^ and Shirley A. Campbell^ 

SUMMARY 

A line-by-line radiative transfer program that simulates 
the response of a Gas Filter Correlation Radiometer (GFCR) as 
a function of altitude is presented and discussed. The program 
was developed to specifically solve downward viewing GFCR 
simulation problems, but can be easily adapted to many 
different infrared sensor requirements. The overlaid structure 
and specific task subroutines permit desired modifications 
to be made with ease. 

The program performs the task of calculating atmospheric 
transmittance and upwelling radiance as a function of wave- 
number, altitude, and primary gas concentration, and writes 
the results on a temporary storage disk. The program then uses 
these results to calculate the output response of the GFCR. 

Both input and output of the program are described, as well 
as a sample test case. The FORTRAN variables, sxibroutines , and 
overlaid programs are listed and explained. A flowchart is 
also furnished. 


^ Research Associate, Department of Physics and Geophysical 
Sciences, Old Dominion University, Norfolk, VA 23508. 

^ Research Associate, Department of Physics and Geophysical 
Sciences, Old Dominion University, Norfolk, VA 23508. 



INTRODUCTION 


The use in recent years of, high sensitivity and high 
effective spectral resolution (less than 0.1 cm“M instru- 
mentation in the remote measurement of trace atmospheric gases 
has required refinements in the computational methods utilized 
in evaluating the atmospheric transmittance in the infrared 
spectral region. The principle of detection of trace atmos- 
pheric gases, specifically gaseous pollutants, has been 
described by Ludwig (ref. 1) . The success of this principle 
not only depends upon the ability of the instrument technique 
to discriminate between the pollutant and interfering spectral 
lines, but also requires an accurate knowledge of high resolu- 
tion pollutant and atmospheric spectra. This can Only be 
accomplished by applying a line-by-line atmospheric radiative 
transfer program. 

The high resolution of the GFCR instrument requires the 
usage of a line-by-line program for the purpose of analyzing 
long path trace atmospheric gas measurements . The program 
presented in this paper is a combination of an improved version 
of the line-by-line program described in reference 2 and several 
desirable characteristics from other line-by-line radiative 
transfer programs (refs. 3 and 4) . The objective was to develop 
an efficient, generalized line-by-line atmospheric radiative 
transfer program, which could be readily adapted to many 
different infrared sensor research needs. The program currently 
accomodates nadir viewing sensors but can be modified for limb 
scan sensors via the appropriate geometrical interpretation of 
atmospheric layer transmission values. This paper describes 
the program and the basic analytical concepts upon which it is 
based, as well as its operation. 
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SYMBOLS 

concentration of absorbing gas, parts per million 
(ppm) 

Boltzmann distribution constant, cm Kelvin 

monochromatic upwelling radiance, watts cm“^sr“^ 

energy of the lower state, cm~^ 

Chapman function, dimensionless 

sensor altitude index, dimensionless 

uppermost layer altitude index, dimensionless 

waveniamber dependent sun irradiance at the top 
of the atmosphere, watts cm“^sr“^cm“^ 

atmospheric layer index, dimensionless 

thickness of layer, cm 

constant, 1.5 for water vapor and 1.0 for all 
other molecules, dimensionless 

Planck blackbody radiation, watts cm'^sr"^ 

mean pressure of layer, atm (1 atm = 1.01325 E+5 
Newtons/me ter^ ) 

equivalent pressure, atm 

adjusted spectral line intensity, atm"^cm~^ 
spectral line intensity, atm“^cm“^ 
layer temperature, Kelvin 

reference temperature corresponding to spectral 
line parameters, Kelvin 

surface temperature , Kelvin 

altitude index, dimensionless 

adjusted spectral line half-width, cm“^ 

line half-width, cm~^ 
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6. Lorentz line shape, dimensionless 

^in 

e wavenumber dependent surface emittance, dimensionless 

<. wavenumber dependent absorption coefficient, 

atm"^cm“^ 

T gaseous transmittance at a particular altitude, 

dimensionless 

0 sun zenith angle, degrees 

w wavenumber (inverse wavelength) , cm“^ 

w averaged waveniamber, cm“^ 

SUBSCRIPTS : 

1 gas species number 

n spectral line number 

DICTIONARY OF FORTRAN VARIABLES 


ABSCOFl 

Absorption coefficient for primary gas at 
OMEGA ( (cm“^atm“M /cm” M • 

ABSC0F2 

Absorption coefficient for interfering gases 
at OMEGA ( (cm” ^atm'M /cm“^ ) . 

ABSCOFT 

Total absorption coefficient for all active 
infrared gases (ABSCOFl + ABSC0F2) ( (cm”^ 
atm" ^ ) /cm”^ ) . 

ADJALPH 

Adjusted half-width of WZ (cm”^). 

ALPHA 

Half-width value at half maximum of WX from 
spectral line parameter tape (cm“M • 

ALPHAD (I,NLAY) 

Half-width values at half maximum of WZ for 
a Doppler spectral line for species I (cm"^) 

ALPHAL (I,NLAY) 

Adjusted half-width values at half maximum 
of WZ for a Lorentz spectral line for 
species I (cm”^). 

ALPHAV (I,NLAY) 

Adjusted half-width values at half maximum 
of WZ for a Voigt spectral line for species 
I ( cm " M . 
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ATMTAU 

(NLAY ,NXM) 

For atmospheric layers, calculated trans- 
mission of atmosphere at the top of the layer 
at OMEGA; for instrumentation and calibration 
layers, calculated transmission of that layer 
at OMEGA. 

BBCNEW 

(NLAYI ) 

Radiation of the cold blackbody calibration 
source attenuated by the vacuum cell at OMEGA. 

BBCNEWC 

(NLAYI) 

Radiation of the cold blackbody calibration 
source attenuated by the calibration cell and 
the gas cell at OMEGA. 

BBCNEWG 

(NLAYI ) 

Radiation of the cold blackbody calibration 
source attenuated by the gas cell at OMEGA. 

BBCNEWV 

(NLAYI) 

Radiation of the cold blackbody calibration 
source attenuated by the calibration cell 
cind the vacuum cell at OMEGA. 

BBCOLD 

(NLAYI) 

Radiation of the cold blackbody calibration 
source attenuated by the vacuum cell at the 
previous OMEGA. 

BBCOLDl 

(NLAYI) 

Temperature of the cold blackbody calibra- 
tion source (Kelvin) . 

BBCOLDC 

(NLAYI) 

Radiation of the cold blackbody calibration 
source attenuated by the calibration cell 
and the gas cell at the previous OMEGA. 

BBCOLDG 

(NLAYI) 

Radiation of the cold blackbody calibration 
source attenuated by the gas cell at the 
previous OMEGA. 

BBCOLDV 

(NLAYI) 

Radiation of the cold blackbody calibration 
source attenuated by the calibration cell and 
the vacuum cell at the previous OMEGA. 

BBHNEW 

(NLAYI ) 

Radiation of the hot blackbody calibration 
source attenuated by the vacuum cell at OMEGA. 

BBHNEWC 

(NLAYI) 

Radiation of the hot blackbody calibration 
source attenuated by the calibration cell 
and the gas cell at OMEGA. 

BBHNEWG 

(NLAYI) 

Radiation of the hot blackbody calibration 
source attenuated by the gas cell at OMEGA. 

BBHNEWV 

(NLAYI) 

Radiation of the hot blackbody calibration 
source attenuated by the calibration cell 
cUid the vacuum cell at OMEGA. 
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BEHOLD (NLAYI) 

BBHOLDC (NLAYI) 

BBHOLDG (NLAYI) 

BBHOLDV (NLAYI) 

BBHOT (NLAYI) 

BLCKIN (NLAYI) 

BRDFAC 

BROAD 

C 

C7 

C8 

C9 

CHAP ZEN 
CMINV 

COMPABS (I,NLAY) 

COSS 

DELTAW 

DIST 


Radiation of the hot blackbody calibration 
source attenuated by the vacuum cell at the 
previous OMEGA. 

Radiation of the hot blackbody calibration 
source attenuated by the calibration cell and 
the gas cell at the previous OMEGA. 

Radiation of the hot blackbody calibration 
source attenuated by the gas cell at the 
previous OMEGA. 

Radiation of the hot blackbody calibration 
source attenuated by the calibration cell 
and the vacuum call at the previous OMEGA. 

Temperature of the hot blackbody calibration 
source (Kelvin) . 

Instrument's internal blackbody temperature 
(Kelvin) . 

Primary gas pressure broadening factor. 

Primary gas pressure broadening coefficient. 

2 * CMINV. 

First radiation constant (1.1908E-12 watts 
cm~^ sr cm"^ ) . 

Boltzmann's constant (1.439 cm K) . 

Doppler coefficient (6.7675E-8). 

Exponential component of the Chapman 
transmission function of the atmosphere for 
the solar zenith angle . 

Distance above and below OMEGA considered 
as center line absorption in the transmittance 
calculation at wavenumber OMEGA (cm“^ ) . 

Line profile component of absorption 
coefficient at OMEGA. 

Cosine of the solar zenith angle. 

Minimum allowable integrating waveniimber 
increment (cm“^ ) . 

Width of the spectral band pass. 
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DV ( JJ ,NEMIS ,NSURF) 

Differential signal resulting from the 
transmission of extemal energy between the 
gas and vacuum cell (watts cm"^ sr~^ ) . 

EL 


Energy of lower state of transition of WZ 
from spectral line parameter tape (cm“^ ) . 

EMISBB 


Emissivity of the instr\iment internal 
blackbody . 

EMISS (NEMIS) 

Thermal emissivity of the earth's surface 
within the effective field of view. 

ERFC 


Error function describing the correction 
for the Voigt profile approximation. 

ERRDN (NLAYI) 

Error in the signal resulting from an 
imbalance between the optical paths (%) . 

FCD 


Doppler intensity function. 

FCERR 


Error function for the Voigt intensity 
approximation . 

FCL 


Lorentz intensity function. 

FCV 


Voigt intensity function. 

FILTERW 

(NOG) 

Instrument filter wavenumber at which input 
transmission values are located (cm“^). 

GAS ANEW 

(JJ, NEMIS, 
NSURF) 

Total upwelling atmospheric radiation 
attenuated by the gas cell at OMEGA. 

GASAOLD 

(JJ, NEMIS, 
NSURF) 

Total upwelling atmospheric radiation 
attenuated by the gas cell at the previous 
OMEGA. 

GASBNEW 

(NLAYI) 

Radiation of the internal blackbody 
attenuated by the gas cell at the previous 
OMEGA. 

GASBOLD 

(NLAYI) 

Radiation of the internal blackbody 
attenuated by the gas cell at the previous 
OMEGA. 

GASCONC 

(iline{l}, 

NLAY) 

Concentration of gas ILINE (L) in layer 
(ppm) . 

GASENEW 

(JJ, NEMIS, 
NSURF) 

Total upwelling atmospheric radiation 
attenuated by the gas cell including gas 
cell thermal emission at OMEGA. 
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GASEOLD (JJ,NEMIS, 
NSURF) 

GBALDIF (NLAYI) 

GCLTAIN (NLAYI) 
IDENT (16) 


ILINE (NOMEG) 

IPOINT 

IPOLLUT 

IS 

JCOUNT 
JPROF (NLAY) 

KBI 

MOLWT(IS) 

MOLWTIS 

NEMIS 

NLAY 


Total upwelling atmospheric radiation 
attenuated by the gas cell including gas 
cell thermal emission at the previous 
OMEGA. 

Difference of the radiation from the hot 
and cold blackbody balance sources 
attenuated by the gas cell path. 

Average transmission of the gas cell for 
the internal blackbody radiation. 

Identification numbers corresponding to 
active infrared gas. 

1 - Carbon Monoxide (CO) 

2 - Water (H 2 O) 

3 - Sulfur Dioxide (SO 2 ) 

4 - Ammonia (NH 3 ) 

5 - Methane (CH 4 ) 

7 - Nitrous Oxide (N 2 O) 

10 - Nitric Oxide (NO) 

12 - Carbon Dioxide (CO 2 ) 

Array of identification nvimbers corres- 
ponding to spectral lines in OMEGSTR. 

Working pointer for OMEGSTR. 

Identification number of primary gas . 

Identification n-umber of species WZ from 
spectral line parameter tape. 

Number of integration steps. 

Line profile array: 1 if Lorentzian, 

2 if Voigt. 

Number of atmospheric layers. 

Array of molecular weights corresponding 
to the IDENT array. 

Molecular weight of species IS (grams/ 
molecules) . 

Niimber of earth's emissivities to be 
considered. 

Total number of input homogeneous layers 
(atmospheric + instrument + calibration) . 
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NLAYI 
NOC 
NOMEG 
NS PEC 

NSTORW 
NS UN 
NSURF 


Number of instrument layers. 

Number of input FILTERW values. 

Nvunber of spectral lines in OMEGSTR. 

Number of active infrared gases being 
considered. 

Number of spectral lines in WSTOR. 

Number of input SUNW values. 

Niomber of surface temperatures to be 
considered. 


NXM 

OMEGA 


OMEGSTR (NOMEG) 


Nvimber of primary gas concentration 
multipliers (XMULT) . 

The center waven limber of the sub interval 
being considered for transmittance 
calculations (cm~i ) . 

Array of spectral lines from OMEGA-CMINV 
to OMEGA + CMINV. 


OPATH (NLAY) 
PLANCK 

PLANK (NLAYI) 

PRES (NLAY) 

PROF AC 1 (I, NLAY) 
PROFAC2 (I, NLAY) 

RADATM (NLAY, NXM) 


Optical path of layer (atm-cm) . 

Thermal emission of layer at OMEGA (watts 

cm~^ sr “M . 

Thermal emission of the instrument's 
internal blackbody at OMEGA (watts 
cm“^ sr“ ^ ) . 

Pressure of layer (atm) . 

Multiplicative variable for line profile. 

Multiplicative variable for line profile 
( I W-Wo 1 /ALPHA) . 

For atmospheric layers, the radiation 
upwelling from the atmosphere as a result 
of atmospheric gaseous molecular emission 
at OMEGA; for instrument and calibration 
layers, the radiation as a result of 
gaseous molecular emission of that layer 
at OMEGA ((watts cm“^sr "^ ) /cm"M • 
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RADCOLC 

(NLAYI) 

Integrated radiation of the cold blackbody 
calibration source attenuated by the 
calibration cell and the gas cell. 

RADCOLD 

(NLAYI) 

Integrated radiation of the cold blackbody 
calibration source attenuated by the 
vacuum cell. 

RADCOLG 

(NLAYI) 

Integrated radiation of the cold blackbody 
calibration source attenuated by the gas 
cell . 

RADCOLV 

(NLAYI) 

Integrated radiation of the cold blackbody 
calibration source attenuated by the 
calibration and vacuum cells. 

RADGASA 

(JJ,NEMIS, 

NSURF) 

Integrated total upwelling atmospheric 
radiation attenuated by the gas cell. 

RADGASB 

(NLAYI) 

Integrated radiation of the internal 
blackbody attenuated by the gas cell. 

RADGASE 

(JJ,NEMIS, 

NSURF) 

Integrated total upwelling atmospheric 
radiation attenuated by the gas cell 
including gas cell thermal emission. 

RADHOT 

(NLAYI) 

Integrated radiation of the hot blackbody 
calibration source attenuated by the 
vacuum cell. 

RADHOTC 

(NLAYI) 

Integrated radiation of the hot blackbody 
calibration source attenuated by the 
calibration cell and the gas cell. 

RADHOTG 

(NLAYI) 

Integrated radiation of the hot blackbody 
calibration source attenuated by the gas 
cell . 

RADHOTV 

(NLAYI) 

Integrated radiation of the hot blackbody 
calibration source attenuated by the 
calibration cell and the vacuum cell. 

RADOLD 1 

(NLAY,NXM) 

For atmospheric layers, the radiation 
upwelling from the atmosphere as a result 
of atmospheric gaseous molecular emission 
at the previous OMEGA; for instrument and 
calibration layers, the radiation as a 
result of gaseous molecular emission of 
that layer at the previous OMEGA. 

RADOTOT 

(JJ,NEMIS, 

NSURF) 

Total upwelling radiation at the previous 
OMEGA. 
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II 

RADSOLD (JJ,NEMIS, 
NSURF) 

RADSUN (NLAY,NEMIS, 
NXM) 

RADSUNO (NLAY,NEMIS, 
NXM) 

RADSURP (JJ,NEMIS, 
NSURF) 

RADTOT (JJ,NEMIS 
NSURF) 

RADVACA (JJ,NEMIS, 
NSURF) 

RADVACB (NLAYI) 

RATLV (I,NLAY) 
REFTEMP (I) 

RESPG (NOG) 

RESPIN 
RESPV (NOG) 

ROOTEMP (NLAY) 

S 

SIGAV (JJ,NEMIS, 
NSURF) 

SIGBV (NLAYI) 


Radiation upwelling from the atmosphere 
as a result of earth surface emission at 
the previous OMEGA. 

Radiation upwelling from the atmosphere as 
a result of incident solar energy at 
OMEGA ((watts cm“^ sr ~^ ) /cm ~^ ) . 

Radiation upwelling from the atmosphere 
as a result of incident solar energy at 
the previous OMEGA. 

Radiation upwelling from the atmosphere 
as a result of earth surface emission at 
OMEGA ((watts cm'^sr /cm "^ ) . 

Total upwelling radiation at OMEGA 
(RADSUN + RADSURF + RADATM) ( (watts 
cm“2sr /cm~M • 


Integrated total upwelling atmospheric 
radiation attenuated by the vacuum cell. 

Integrated radiation of the internal 
blackbody attenuated by the vacuum cell. 

Ratio of Lorentz to Voigt half-widths 
for species I. 

Reference temperature corresponding to 
identification numbers of spectral line 
parameters from tape (Kelvin) . 

Normalized transmission values for gas 
cell given at FILTERW. 

Instrument detector response. 

Normalized transmission values for 
vacu^am cell given at FILTERW. 

Square root of layer temperature. 

Adjusted line strength of WZ (atm”^cm“^). 

Signal resulting from external blackbody 
radiation attenuated by the vacuum cell 
and the aperture (watts cm~^sr“M . 

Signal resulting from the internal 
blackbody radiation attenuated by the 
vacuum cell and the aperture (watts 
cm“^ sr“^ ) . 
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SUNCOM (NEMIS) 
SUNFLUX (NSUN) 


SUNINTP 
SUNW (NSUN) 

SURFCOM (NEMIS, 
NSURF) 


Component of solar energy unattenuated 
by the atmosphere (watts cm~^sr“^). 

Wavenumber dependent solar energy incident 
at the top of the atmosphere (watts 
cm~^sr ) . 

Interpolated solar energy at OMEGA. 

Waveniambers corresponding to the SUNFLUX 
input values. 

Component of surface energy contribution 
unattenuated by the atmosphere (watts 
cm~^sr ) . 


SZ 


Line strength values of WZ from spectral 
line parameter tape (atm~^cm~^) . 


TAU (NXM) 
TAUA (NLAYI) 


TAUGFIL 


Transmission of layer at OMEGA. 

Calculated average transmission coefficient 
for the vacuum cell aperture required to 
oalance the optical paths for two external 
blackbody temperatures. 

Interpolated transmittance of gas cell 
filter at OMEGA. 


TAUOLD ( NLAY , NXM) 


TAUSLAT (NXM) 
TAUVFIL 


For atmospheric layers, the calculated 
transmission of the atmosphere at the 
top of the layer at the previous OMEGA; 
for instriament and calibration layers, 
the transmission of that layer at the 
previous OMEGA. 

Total atmospheric slant path coefficient. 

Interpolated transmittance of vacuum 
cell filter at OMEGA. 


TCONSQ 


TEMP (NLAY) 
TEMPCON 


TERFC 


Temperature dependence of the rotational 
partition function; (TEMPCON) ^ unless 
water, which is (TEMPCON) . 

Temperature of layer (Kelvin) . 

Reference temperature of gas/temperature 
of layer (temperature adjustment for line 
half -width) . 

Error function term for Voigt profile. 
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TERFC2 
THICK (NLAY) 

TINV (NLAY) 

TOTATM (NLAY,NXM) 


TOTRAD (JJ,NEMIS, 
NSURF) 

TOTSUN ( NLAY , NEMIS , 
NXM) 

TOTSURF (JJ, NEMIS, 
NSURF) 


TRANS (NLAY, NXM) 


TSURF (NSURF) 


V (JJ, NEMIS, NSURF) 


VACANEW ( J J , NEMIS , 
NSURF) 

VACAOLD (JJ, NEMIS, 
NSURF) 


VACBNEW (NLAY I) 


VACBOLD (NLAYI) 


TERFC squared. 

Thickness of layer (cm) . 

Inverse of layer temperature. 

For atmospheric layers, the integrated 
radiation upwelling from the atmosphere 
as a result of atmosphere gaseous molecular 
emission; for instrument and calibration 
layers, the integrated radiation as a 
result of gaseous molecular emission of 
that layer . 

Total integrated upwelling radiation. 


Integrated radiation upwelling from the 
atmosphere as a result of incident solar 
energy. 

Integrated radiation upwelling from the 
atmosphere as a result of earth surface 
emission . 

For atmospheric layers, the calculated 
integrated transmission of the atmosphere 
at the top of the layer; for instrument 
and calibration layers , the integrated 
transmission of that layer. 

Temperature of the earth's surface 

within the effective field of view (Kelvin) . 

Average signal resulting from the transmission 
of external energy between the gas and 
vacuum cells (watts cm~^sr~^). 

Total upwelling atmospheric radiation 
attenuated by the vacuum cell at OMEGA. 

Total upwelling atmospheric radiation 
attenuated by the vacuum cell at the 
previous OMEGA. 

Radiation of the internal blackbody 
attenuated by the vacuum cell at OMEGA 
(watts cm~^sr~M . 

Radiation of the internal blackbody 
attenuated by the vacuum cell at the 
previous OMEGA (watts cm~^sr”M • 
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VBALDIF (NLAYI) 


VGTPAR 

VGTPAR2 

VGTPAR4 

W3 

WI 

WF 

WRATLV 


Difference of the radiation from the hot 
and cold blackbody balance sources 
attenuated by the vacuum cell path. 

A Voigt parameter. 

Reciprocal of the Voigt parameter squared. 
VGTPAR squared. 

OMEGA cubed. 

Initial wavemamber used for transmittance 
calculation (cm"M • 

Final wavemamber used for transmittance 
calculation (cm“M . 

Ratio of Lorentz to Voigt half-widths 
for WZ. 


WRATVl 


1 - WRATLV. 


WSTOR (NSTORW) 


WZ 

XMULT (NXM) 


Array of spectral lines, the first of 
which is £ OMEGA and the remainder are ^ 
OMEGA but £ OMEGA + CMINV. 

Wavenumber values from spectral line 
parameter tape (cm“M . 

Primary gas concentration multipliers. 


ZENITH 


Sun zenith angle (degrees) . 


PROBLEM DEFINITION AND METHOD OF SOLUTION 


A line-by-line radiative transfer computer program was 
needed to efficiently perform the task of data reduction for 
several versions of the GFCR which are being developed by 
NASA/Langley Research Center (LaRC) . The writing of the SMART 
program was undertaken to minimize data reduction cost by the 
efficient utilization of computer storage through overlaying, by 
reducing execution time, and by combining several phases of the 
previous data reduction procedure, i.e., instrument balance and 
calibration (ref. 1) and data signal simulation. The modular 
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structure of the program was designed to permit modifications of 
the computational algorithms without affecting the program frame- 
work. The individual overlaid programs and subroutines are 
explained in the Program Organization and Description section. 

In addition, thorough documentation was necessary for program 
clarity. 

In order to minimize computation time without a loss in 
accuracy, certain assumptions were made concerning the atmos- 
phere and radiative transfer processes in the algorithms used 
to perform transmittance calculations. To compute the trans- 
mittance from the earth's surface to a sensor altitude, h , 
the modeled atmosphere between these points is divided into 
a number of homogeneous layers, i.e., regions within which 
the temperature, total pressure, and concentrations of the 
primary and interfering molecular absorbers are uniform. The 
error in this approximation may be made as small as desired 
by subdividing the atmosphere into a sufficiently large number 
of layers. The total radiance incident on a sensor at altitude 
h is given by 

E(h) = E(to)do3 ( 1 ) 

where Ato is the spectral bandpass of interest and E(uj) is 
the total monochromatic upwelling radiance. 

For a cloud-free, non-scattering atmosphere under local 
thermodynamic equilibrium, the atmospheric radiative trans- 
fer equation for the total monochromatic upwelling radiant 
energy, E(to) , as viewed by a nadir type of sensor can be 
written as 
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E(w) = e(o)) N°(o),T ) x(w,h) 

s 

^ N°(.,T(Z)) d. 

( 2 ) 

+ i [1 - e(o))] cos 9 Hg(oj) 

• [T(0),h)] [T(W,h')] f(6) 


where e(o)) is the wavenumber dependent surface emittance, 

N°(6o,T ) is the Planck blackbody function which is dependent on 
s 

wavenumber and surface temperature, T , or radiating gas 

W 

temperature at a particular altitude, T(z) . The monochromatic 
transmittance of the atmosphere between the emitting surface z 
and the altitude of the sensor, h , is represented by T(w,h) , 
and the monochromatic vertical transmission of the entire 
modeled atmosphere is represented by T(o),h') . The solar zenith 

angle is 0 and the wavenumber dependent sun irradiance at the 
top of the atmosphere is H . The Chapman function (ref. 1) , 
f(0) , is equal to sec 9 for 0° £ 9 £ 60° and is equal to 
the Chapman polynomial for 9 > 60° . The three terms on the 
right-hand side of equation (2) represent, respectively, the 
earth's surface emission, the atmospheric emission, and the 
solar reflected energy (see figure 1) . All of these components 
must be considered in the solar-thermal overlap region at 
4.6 ym. These terms are represented in SMART by RADSURF, 

RADATM, and RADSUN, respectively. 

A closer examination of the form of equations (1) and (2) 
indicates two very distinct differences between the methodology 
of SMART and other conventional line-by-line radiative transfer 
programs. The first difference is the order in which the 
integrals over Au and over the change in altitude. Ah , are 
performed. By initially integrating over Ah monochromatically , 
the inflexibility of transmittance averaging is eliminated, i.e.. 
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Figure 1 


Interaction of radiation with the atmosphere. 


the monochromatic transmission at any altitude h is given 

by 


h 

T (( 0 ,h) = n X (O) 
z=l 

and not by 

?(Ao,h) = ^ E 

Ao) 

The spectral resolution of the transmittance values, i.e., the 
number of w's at which the absorption coefficient is calculated, 
is defined by the user in subroutine GETW. The algorithm used 
for our CO calculations is explained later in this section. 

The second difference between SMART and most other line-by- 
line radiative transfer programs is the specification of the order 
of limits of integration over altitude. The lower limit of inte- 
gration is the radiating source and the upper limit is the sensor 
altitude which allows for practical computation of the mono- 
chromatic and total integrated transmission at the top of each 
atmospheric layer as seen by equation (3) , This eliminates 
redundant calculations of atmospheric transmittance in evaluating 
signals from aircraft platform sensors at various altitudes. 

Theoretically, the total absorption coefficient at any 
wavenumber, m , consists of contributions from all spectral 
lines; however, for practical computational purposes, only 
lines within the vicinity of m are considered for calculation 
in SMART. The contribution to the absorption coefficient of 
lines in the vicinity of u can be divided into two parts, 
direct and wing. Those lines which lie within an interval 
defined by approximately 100 half-widths of the primary gas to 
each side of m are considered to be the direct contributors 


z) 


n T (o), z) 

Z = 1 


(3) 


( 4 ) 
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to the absorption coefficient, while those lines lying outside 
of this interval result in wing absorption. In the SMART 
calculation of direct line absorption contribution, at each 
incremented step a symmetric interval (CMINV) is considered 
about the center wavenumber. This interval is constant through- 
out the entire band. For atmospheric carbon monoxide using the 
Lorentzian line profile, an interval of 5 cm~^ was selected since 
this value was approximately 100 times the line half-width of 
CO. For spectral lines beyond 5 cm“^, the direct contribution 
to the absorption coefficient is very small (ref. 8). Presently, 
the SMART program does not include wing contribution. 

The SMART program is divided into two primary level over- 
lays (see figure 2) . The first lower level overlay, as described 
in the Program Organization and Description section, performs 
initialization of all variables required for the calculation of 
the monochromatic gas transmission (TAU) for each atmospheric, 
instrument cell, and instrument calibration layer. The procedure 
for the computation of TAU is initiated by a call of the READTP 
subroutine, which reads all spectral reference information; i.e., 
line position (WZ) , half-width (ALPHA) , intensity (SZ) , lower 
energy level (EL), and species identification n\imber (IS) from 
the spectral line parameter reference tape for the band interval 
WI - CMINV to WF + CMINV. To analyze the GFCR sensor carbon 
monoxide data measurements in the 4.6 ym region, the McClatchey 
spectral line parameter tape (ref. 5) was used to obtain all 
spectral reference information. The parameters read should be 
under room temperature and standard pressure (296.0 K and 1 atm) 
conditions or the appropriate conversion must be performed prior 
to reading these parameters in subroutine READTP (see Appendix C) . 

The gaseous transmittance at a particular altitude as a 
function of wavenumber is given by 


TAU(io) = t(w) = exp ^ Kj^(u))p 


(5) 
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Figure 2 


SMART program structure 













where k^(u) is the waveniomber dependent absorption coefficent 
of gas species i (their sum being ABSCOFT) , p (PRES) is 
the total mean pressure of layer k , (GASCONC) is the con- 

centration of absorbing gas i , and % (THICK) is the thick- 
ness of the kth layer. 

In general, the absorption coefficient of the nth line 
of species i is described by 


■'in'”* = ®inSin<“> 


( 6 ) 


where is the layer temperature (TEMP) corrected line 

intensity and is the line shape function, i.e., Lorentz, 

Voigt, or Doppler. The Lorentz line shape is given by 


^in(“> = ¥ 


“in 


( 0 ) - 0 ).^) 




(7) 


where (ADJALPH) is the temperature and pressure dependent 

line half-width for layer k and is calculated by 

/T \ V2 

ADJALPH = = aoPeV'T"/ * 

The reference line half-width, (ALPHA) , is read from the 

spectral line parameter tape at temperature T^ (REFTEMP) . 

The equivalent pressure, p , is a function of the ratio of 
self-broadening to the nitrogen-broadening efficiency (BROAD) , 
the total pressure (PRES) , and the concentration of the 
absorbing gas (GASCONC) as given by 

Pg = [GASCONC * (BROAD - 1) +1] *PRES . (9) 

In the case of trace gases in the atmosphere, the equivalent 
pressure is set equal to the total atmospheric pressure since 
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self-broadening is insignificant. For instrument and calibra- 
tion layers, self-broadening is assumed and the corresponding 
equivalent pressure for BROAD is calculatedi 

The line intensity (S) depends upon the temperature 
through the Boltzman distribution factor (ref. 6) and is 
expressed as 


S . = 
in 





( 10 ) 


or 


. S = SZ *TC0NSQ * EXP [-C8* (TEMPCON - 1) *EL/REFTEMP] (11) 

where m is 1.5 for water vapor and ozone and 1.0 for other in- 
frared active molecules, such as CO, CO 2 , N 2 C, and C8 is the Boltzman 
distribution constant. At this point, the user has the option of 
calculating ABSCOFT as determined by the Lorentz function or as 
determined by an approximation to the Voigt function detailed by 
Kielkopf (ref. 7). The flexibility of SMART allows easy insertion 
of the Doppler line profile calculation by the addition of an 
option as defined by JPROF in the same program location as the 
Lorentz and Voigt function. The line shape for each layer (JPROF) 
must be designated in the input and need not be the same for each 
layer. The direct contribution to the absorption coefficient from 
all infrared active lines whose center positions lies within CMINV 
of the wavenumber under consideration (OMEGA), i.e., the • interval 
OMEGA - CMINV to OMEGA + CMINV, constitutes the total direct 
absorption at OMEGA. For additional spectral line information, 
the interval of direct contribution can easily be extended by 
increasing the dimensions of the appropriate variables. The 
total absorption at OMEGA for this interval is determined by 
the summation of the absorption coefficients for all lines 
within the interval, in addition, the user can automatically 
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I 

obtain the total absorption coefficient (ABCOFT) that 
corresponds to a maximum of ten different primary gas vertical 
mixing ratio profiles as defined by XMULT. Each value of 
XMULT is multiplied by the volume mixing ratio of the primary 
gas in every layer resulting in a bias shifting of the input 
vertical mixing ratio profile. ABSCOFT is calculated as the 
sum of two components, ABSCOFl, which is the absorption 
coefficient for the primary gas, and ABSC0F2, which is the 
absorption coefficient for all other interfering gases. The 
absorption effects of continue, such as nitrogen and water 
vapor, could easily be considered in the form of an added 
third term to ABSCOFT after the necessary algorithms for 
accurately calculating the continuum absorption in the spectral 
region under consideration have been developed. 

After completion of the monochromatic calculation of ABSCOFT 
at OMEGA for the first layer, the transmittance is calculated as 
shown by equation (5). This procedure is repeated resulting in 
a TAU value for each layer k . The self-emission (RADATM) of 
each layer is then calculated via the wavenumber dependent energy 
described by the Planck blackbody function used in subroutine 
PLANKF. The temperature used in PLANKF is the mean temperature 
(TEMP) of the emitting layer. The emissivity (EMISS) of the 
layer is determined by the calculated transmittance of that layer 
as given by 

EMISS = 1 - TAU . (12) 

The total monochromatic transmittance (ATMTAU) at the top 
of each atmospheric layer is then calculated as 

ATMTAU(o)) = TAUj^(w) (13) 

and written on temporary storage disk 6 with the corresponding 
wavenumber and self-emission component. 
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The next operation performed by SMART is the determination 
of the wavenumber increment step size for the pvirpose of 
repeating the transmittance calculations for the entire spectral 
band being considered, i.e., WI - CMINV to WF + CMINV. This 
task is accomplished by the GETW subroutine. The size of the 
wavenumber increment used in absorption coefficient calculations 
is arbitrary and depends upon the spectral resolution required 
by the user, since the absorption coefficient varies rapidly as 
a function of relative position to the spectral line centers. 
Subroutine GETW can be replaced by any of several schemes (refs. 
4, 8, and 9) proposed for selection of wavenumber positions for 
line-by-line absorption calculations. For applications of 
SMART to the GFCR carbon monoxide sensor, a minimum mesh size 
(DWMIN) of .01 cm“^ and a maximum mesh size of .5 cm~\ which is 
1/8 of the average distance between CO spectral lines, was 
determined to most efficiently describe the wavenumber dependent 
absorption coefficient. When the spectral line density being 
examined is high, a large number of integration points are 
calculated; however, if the spectral line density is low, fewer 
points of integration are considered. The transmittance cal- 
culation for a new OMEGA is performed and the corresponding 
results, i.e., ATMTAU (w) and RADATM (u) , are written on 
temporary storage disk 6. This procedure is repeated until 
ATMTAU and RADATM for the entire spectral band are written. 

The second lower level overlay is comprised of four 
sections. The first section of the overlay reads the input 
data required for computation of instrument parameters, and 
the remaining variables of the RADSURF and RADSUN contribu- 
tions given by equation (2) . After initialization of variables, 
subroutine FILTER? is called to evaluate by interpolation the 
wavenumber dependent variables, i.e., instrument filter 
functions and solar irradiance. values. The second section 
calculates the earth's monochromatic emission (RADSURF) as a 
function of one or more surface temperatures and emissivities . 
This task is accomplished by employing subroutine PLANKF in 
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the calculation of the wavenumber dependent Planck distribution 
of energy. 

The solar contribution (RADSUN) is calculated in accordance 
with equation (2) by subroutine SUNNY/ where the Chapman 
function is employed in determining the slant path transmittance 
of incident solar radiance. The three terms of equation (2) , 
i.e./ RADSURF, RADATM, and RADSUN/ are then obtained by attenu- 
ating each component by the appropriate atmospheric trans- 
mittance (ATMTAU) value for an altitude corresponding to the 
top of each atmospheric layer. 

The third section calculates the internal radiance values 
for the GFCR sensor and the necessary parameters for balance and 
calibration. [A detailed description of instruments based on 
the gas filter correlation technique is presented in references 
1 and 2] . In addition, the simulated signals (DV and V) 
corresponding to the sensor detection of total upwelling radiance 
(RADTOT) are calculated. For convenience, the trapezoidal rule 
is used for integration of all variables over the spectral 
band being examined. This process is performed by the subroutines 
XNTEG, XNTEG2, XNTEG3, and XNTEG4 . The dimensions of the 
variables to be integrated determines which of the four sub- 
routines is called. The required output parameters, both 
atmospheric and sensor, are then stored and are readily 
accessible for output listing in section four. 

PROGRAM ORGANIZATION AND DESCRIPTION 


The SMART program is written in FORTRAN IV for the Control 
Data Cyber series computer systems and is overlaid. For our 
sample test case using 17 input layers, the storage required 
is 130300 octal words. The storage is dependent upon the 
number of atmospheric and instrument layers required for 
calculations. Each additional layer requires approximately 
2000 octal words. 


25 



Each overlaid program and its subroutines are listed on 
the following pages. For each program or subroutine, the 
required niomber of octal words of storage and an explanation 
of the function are presented. A dictionary of FORTRAN 
variables used was given in a previous section. A flow chart 
showing the logic flow and the interrelation of the various 
programs and subroutines is shown in figure 3 . 
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DO 170 K»l, NIAY 



Figure 3. Flowchart of SMART program. 
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SUBROUTXKE G£T(D 


SUBROUTINE READTP 




Figure 3. Flowchart of SMART program (continued). 
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PROGRAM SRIGTCA 



Figure 3. Flowchart of SMART program (continued) 
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PROGRAM SMART 


SMART, the program executive, has two main functions. 
First, it supervises the execution of the two lower level 
overlays. Second, storage for each array and program constant, 
which is common to each overlay, is set aside. The storage 
required, i.e., program length + buffer length + labeled 
common length + blank common length, is 13241s . 


I 


5 


10 


15 


0VERL*Y(HAIN,0>0I 

PROGRAH SMART Cl NPUT,0UTPUT«TAPE7>TAPE2,TAPE6I 
C *** THIS PROGRAM PERFORMS LINE-BY-LINE RADIATIVE TRANSFER CALCULATIONS FOR 
r NON-HOMOGENEOUS ATMOSPHERIC MODELS AS A FUNCTION OF ALTITUDE 

C *** THE SECOND SECTION SIMULATES THE RESPONSE OF THE GAS FILTER CORRELATION 
C INSTRUMENT AT THE VARIOUS ALTITUDES 

C0HH0N/AB/TENP(17t,PRES(17>»THICK<17l,0PATHI17)»6ASC0NCI20,17)» 

1 IDENT(20),XNULT(10),JPR0FI17I 

comhon/aa/nxh,nlay,nlayi*noc*kbi,nspec,deltaw>cminv«ipollut, 

1 broad, WI*WF 

CALL OVERLAYCAHMAIN, 1,01 
REWIND 6 

CALL 0VERLAY(AHHAIN,2,0) 

STOP 

END 
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PROGRAM TRANSMT 


TRANSMT is the first overlay subordinate to SMART. 

TRANSMT reads the atmospheric profile information and the 
integration processing parameters from cards. From the spectral 
line parameter tape, wavenumbers, half-widths, intensities, 
energies of the lower state, and identification species num- 
bers are read. For our calculations, McClatchey's spectral 
line parameter tape (ref. 5) was used with appropriate unit 
conversions applied. TRANSMT perfo3nns line-by-line radiative 
transfer calculations and writes on temporary storage disk 6 
center wavenumber, transmittance, and radiance as a function of 
altitude and primary gas concentration. Storage required is 
112274a . 
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5 


10 


15 


ZO 


25 


30 


35 


QVERLAY<nAlN,l,0) 

PROGRAM TRANSMT 

DIMENSION REFTEMP(30)»AtPHAL(300>17)#PR0FACl(300»l7)> 

1 PR0FAC2(300#17)»TAUa0lfC0MPABS(300»17)»TlNV(17), 

1 ROOTEMP ( 17)^ ATMTAU(17f 10) > RAOATM C17# 10) > 

Z AtPHAVC 300f 17),RATLVC300#17)>ALPHAO(300#17) 

REAL MOLWTI20)fMOLWTIS 

COMNON/AB/TEMP<17)»PRESC17)f THICK(17)90PATH(I7),GASCQNC(20,17)# 

1 IDENT(20)#XMULTI10)>JPR0Ff 17) 

C0HH0N/AA/NXM,NLAY»NLAYI,N0C#KBIf NSPEC»DELTAW»CMINV#lPOLLUTf 
1 BROAOfWlfWF 

COMMON WST0RC300) •0NEGSTRI300)>UINE(300) 

DATA C7/1.1906E-12/»C8/1.439/»C9/6.7675E-8/»C10/109./,MSaN/41/ 
DATA (REFTEMP-30*296.) 

0ATA(H0LWT«28.,18.f 64.,17.#16.,4^*,94.#28.,17.»30.»A6.,A4>^^«,6A., 
1 A6.»44.»17.#30.) 

C ♦♦♦INPUT PARAMETERS READ FROM CAROS ♦♦♦ 

READ lOlOf NlAY»NLAYI,NOC>NXM,NSPeC*IOENT 
1010 F0RMATC5I5»5X#20I2) 

READ lOaOt (JPR0F(K|,K«1»NLAY) 

1020 F0RMAT(A0I2) 

READ 1030» WI»WF»DELTAW»CM1NV»1P0LLUT»BR0A0 
1030 FORMATCAFlO.Af I5»F10.4) 

READ lOAOt (XMULT(J)»J-1»NXM) 

10^0 F0RMAT(8F10.4) 

READ 1050» CTEMPC J)f PRES< J)»THICK( J),J«1,NLAY) 

1050 F0RMAT(3F10.4) 

READ 1060# (fGASC0NC(J#K), J-1#16)#K«1,NLAY) 

1060 F0RMATC8E10.3) 

C ♦♦♦ INITIALIZATION ♦♦♦ 

KBI«NLAY-2^NLAYI 

OHEGA-WI 

DWl-0. 

NST0RW«N0MEG»0 

I-O 
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^0 


45 


50 


55 


60 


65 


70 


75 


80 


85 


90 


95 


100 


105 


DO 10 K»1»NIAY 

0PATH(K)«PRES(K)4 THICKIK) 

TINVCK)-1#/TEHPCK) 

ROOTEMP(K)«SORT(TEHP(KM 

BRDFAC»IBR0A0-1« )*6ASC0NCaP0LLUT>K)4l. 

IF (NLAY.EQ.l) PRES < K I -PRES CK ) ♦BRDF AC 
IF (K.GT.KBI) PRES(K)«PRES(KI*BRDFAC 
10 CONTINUE 

C ♦♦♦ READ SPECTRAL LINE PARAMETER TAPE FOR ALL LINES INCLUnEO IN WI-CMINV TO 
C WF ♦CMINV ♦♦♦ 

CALL REA0TP(WI»CMINV»WF»IDENT#NST0RW»N0MEG) 

20 W3-0MEGA ♦OMEGA ♦OMEGA 

C ♦♦♦ DETERMINE SPECTRAL LINES TO BE INCLUDED IN SUBIMTERVAL OF INTEGRATION 
C AND ADJUST CORRESPONDING LINE PARAMETERS ♦*♦ 

IF (OMEGA .EQ. WI ) GO TO 50 

IF tOHEGSTRd) *GE* ( OMEGA-CHINV l> GO TO 45 
DO 25 K«2»N0HEG 
IPOINT-K 

IF (OMEGSTR(K) .GE. (OHEGA*CMINVn GO TO 30 
25 CONTINUE 

30 NOHEG«NOMEG -IPOINT 

I«I -IPOINT +1 
DO 35 J«1»N0MEG 

ILINE(J)«ILINEU ♦IPOINT-II 
OMEGSTRC JI-OMEGSTR( J ♦IPOINT-IJ 
35 CONTINUE 

DO 40 K-lfNLAY 
00 40 J«1#N0ME6 

ALPHALC J»K)»ALPHAL(J ♦IP0INT*1>K) 

PR0FAC1IJ^K)«PR0FAC1(J ♦IP0INT-1#K) 

IF UPROFCKI.EO.ll GO TO 40 
RATLVC J»K|-RATLV( J+IPOINT-lfK) 

40 CONTINUE 

45 IF (WZ «GT. (OMEGA ♦CMINVM GO TO 100 
60 TO 65 
50 I-O 

C ♦♦♦ READ SPECTRAL LINES FOR SUBINTERVAL OF INTEGRATION ♦♦♦ 

REWIND 2 

55 REA0(29l000) WZ » ALPHAS SZ » EL» IS 
1000 FaRMAT(FlO«4>FS*3»E15.0»F12.4»I2l 
IF (eOF(2>) 100#60 

60 IF (WZ *GT« (OMEGA ♦CMINV)) GO TO 100 
65 I»l4l 

IF (OMEGA •EO* WI ) GO TO 70 
NOMEG-NOMEG ♦! 

OMEGSTR(NOHEG)«WZ 
IL!NE(NOMEG)*IS 
NSTORWNSTORW ♦! 

WSTOR(NSTORW) -WZ 

C ♦♦♦ determine LINE PROFILE PARAMETERS AND ADJUST FQR LAYER TERPERATURE ♦♦♦ 
70 00 95 K«lfNLAY 

TEMPCON*REFTEMP( IS) ♦TINV(K) 

TCONSQ-TEMPCON *TEHPCON 

IF (IS .EO. 2) TCONSO- TCONSO ♦SORT ( TEMPCON ) 

AOJALPH-ALPHA ♦PRE$(KI^ SORT ( TEMPCON ) 

S«SZ ♦TCONSO ♦EXP(^C8 ♦(TEMPCON -!•) ♦ EL /REFTEMPI IS ) ) 

IF (JPR0F(K)»E0.1) GO TO 85 
MOLWTIS-MOLWTfISI 

ALPHADn»K)«C9^WZ^R00TEMP(K)^SQRT(28*/M0LWTIS) 

TSTALPH»ADJALPHZALPHAD( I#K) 

ALPHAV(I#K)«.53431^ADJALPH4S0RT(.21687^A0JALPH^ADJALPH 
1 ♦ ALPHAD( I»K)^ALPHAD( IpK)) 

VGTPAR«.83255^ADJALPMZALPHA0(I»K) 

VGTPAR2*1./CVGTPAR^VGTPAR) 

RATLV(I»K)-l./( 1»>A0JALPH^VGTPAR2/ALPHAV( t»Kn 
IF (VGTPAR.lt*!. 5) GO TO 75 
VGTPAR4*VGTPAR2^VGTPAR2 

ERFC*(l.+4.5+V6TPAR2+2.+V6TPAR4l/(1.45.*VGTPAR2 
1 ♦ 3.75*VGTPAR4) 

GO TO 80 

75 TERFC-1./(1.^*47047^VGTPAR) 

TERFC2«TERFC^TERFC 

ERFC«VGTPAR^TERFC4(.61686-.16994^TERFC+1.32554ATERFC2) 
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110 


115 


120 


125 


130 


135 


140 


145 


150 


155 


160 


165 


170 


175 


160 


80 CONTINUE 
65 ALPHUd^K)* AOJAIPH 

IF (JPR0F(K).EQ.2) 60 TO 90 
PROFACKIf K)» ,3183 +S/AOJALPH 
GO TO 95 

90 PROFACl(IfK)*.3183+S*ERFC/ADJALPH 
95 CONTINUE 
GO TO 55 

C ♦♦♦ CALCULATE ABSORPTION COEFFIENT AS DETERMINED BY LINE PROFILE 
100 DO 120 K*1,NLAY 

DO 120 J«1»N0HEG 

PR0FAC2(JfK)-AB$( OMEGA- 0ME6STR U ) ) / ALPHAL ( K ) 

IF f JPPOFfKJ.EO.n GO TO 115 

WRATLV-RATLVtIfKI 

WRATLVl-l.-WRATLV 

FCL-l./ll. + (PR0FAC2U>K»APR0FAC2U#Kn ) 

TSTPROF-ALPHALC JfKi/AtPHAOl J>KI 
IF CTSTPRDF.GT.l) GO TO 105 
EXPTST— .693147*PR0FAC2fJfK)*PR0FAC2( Kl 
IF (EXPTST.lt. -675,84) GO TO 105 
FCO«EXP(EXPTST) 

60 TO 110 
105 FC0*O. 

110 FCERR«( (.8029/PROF AC 2 (J»K) ) - ( . 4207*P ROF AC 2 ( J > K J) I / 

1 (Cl./PR0FAC2UfKn ♦ C PROF AC2 ( J » K ) *PROF AC2 ( J » K )> ♦ 

2(.2030 ♦ U07355*(PROFAC2U#K)*PROFAC2(JtKn ))) 
FCV-WRATLV1*FC0+WRATLV*FCL+ 

1 WRATLV*WRATLV1*FCERR*(FCL - FCO) 

C0HPABSU»K>*6ASC0NC(ILINe(JI,KJ*FCV*PR0FACl(J»K) 

GO TO 120 
115 CONTINUE 

COMPABSU»K)-6ASCONC(ILINE( J)f K)*PR0FAC1IJ»K)/(1.+PR0FAC2( J»K)+ 

1 PR0FAC2(J»K)) 

120 CONTINUE 

DO 125 K*1»NLAY 
DO 125 J«1#NXM 
RADATM(K»J)«0. 

ATMTAUCKf J)«l. 

125 CONTINUE 

DO 170 K»1,NLAY 
IF (K.EO.l) GO TO 135 
DO 130 J-1»NXM 
ATMTAU(K» J)-ATHTAU(K-1» J) 

RADATM(Kf J)*RA0ATH(K-lf J) 

130 CONTINUE 

135 ABSC0F2-ABSC0F1-0. 

IF (K.LE.KBI) GO TO 145 
DO 140 J-1»NXM 
ATMTAU(K,JI-1. 

RADATfKK, JJ-0. 

140 CONTINUE 

145 DO 160 L-1»N0HEG 

IF (ILINE(U .EO. IPOLLUT) GO TO 150 
ABSC0F2- ABSC0F2 • COMPABSCLtK) 

GO TO 155 

150 A8SC0F1-ABSC0F1 ♦COMPASS ( Lf K ) 

155 CONTINUE 

160 CONTINUE 

PLANCK-C7+W3/ (EXP (C8A0MEGA /TEMP (K ) )-l. ) 

ABSCOFT*A0SCOF1 ♦ABSC0F2 

C CALCULATE MONOCHROMATIC ATMOSPHERIC TRANSMITTANCE AND EMISSION AND WRITE 

C RESULTS ON DISC 6 *♦* 

DO 165 J»1»NXM 

IF (K.LE.KBI) ABSCOFT- ABSCOFl *XMULT(J) ♦ABSC0F2 
TAU(J)- EXP(-ABSCOFT AOPATHIK)) 

RADATM(K, J).RAOATM(K> J)*TAU( J)+PLANCK*(1.-TAU(J) » 

ATMTAUlKf J)«ATMTAU(K. J)*TAU( J) 

165 CONTINUE 
170 CONTINUE 

WRITE(6)0MEGA. ( ( ATNTAU( K» J ) > J«l* NXM ) # K-1. NL AY ) » ( ( R AD ATM (K# J ) » 

1 J»l»NXM)#K«lf NLAY)f DWl 

C ♦♦♦ CALCULATE NEW OMEGA ♦♦♦ 

CALL GETW(OMEGAfOELTAWfNSTORW#DWl) 

IF (OMEGA .LE« WF ) GO TO 20 
RETURN 

END 34 



SUBROUTINE GETW 


GETW is a subroutine of TRANSMT. GETW determines the 
waveniimber mesh for integration. The algorithm 1) determines the 
distance from OMEGA to the closest spectral line, 2) cal- 
culates 1/4 of that distance, 3) determines the maximum (DW) 
between DWMIN and the value calculated in step 2, and 4) sets 
the new integration point equal to OMEGA plus DW. The minimum 
integration stepping size is read in TRANSMT. Storage required 
is 1644s . 


1 


5 


10 


15 


SUBROUTINE GETW( W , 0W» NSTORW, DWl ) 

♦♦♦ DETERMINE THE WAVENUMBER MESH FOR INTEGRATION ♦♦♦ 
COMMON WSTOROOO) ,0MEGSTR(300)» ILINEOOO) 
W1«W-WST0R( 1 ) 

WZ-WST0R(2)-U 

DW1-AMAX1(AMIN1(W1>W2>*.25,DWI 
W<W 4DW1 

IF (WST0R(2i .GT. W) RETURN 
K»3 

5 IF (WSTOR(K).GT.Wi GO TO 10 
K«K*1 
GO TO 5 

10 NSTORWNSTORW-K+2 
DO 15 J«1,NST0RW 
15 WST0RUI-WST0R( J + K-2) 

RETURN 

END 
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SUBROUTINE READTP 


READTP is a subroutine of TRANSMT. READTP reads the 
spectral line parameter tape for all lines in the interval 
W1 - CMINV to WF + CMINV and writes them on temporary storage 
disk 2. Storage required is 2000 e . 
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SUBROUTINE READTP C Wl »C H, WF» I DENTt K, L I 

C ♦♦♦ READ SPECTRAL LINE PARAMETER TAPE FOR ALL LINES FROM WI-CMINV TO WF ♦ 

C CMINV 

C VARIABLES INCLUDE LINE LOCATION ( WZ ) ^ HALF-WIDTH (BUfLINE STRENGTH iS), 

C GROUND STATE ENERGY tE)»SPEClES IDENTIFICATION NUMBER (IS) ♦♦♦ 

DIMENSION IDENT(20) 

COMMON WST0Rn00)F0MEGSTR(300l>ILINEC300) 

REWIND 2 

K*1 

L-0 

5 READITflOOO) WZfBLfSfE»IS 

IF <WZ .GT. (WF +CM)> GO TO 25 

IF IWZ .LT, CWI- CM) J GO TO 5 

DO 10 J« 1^16 

IF tIDENTU) .EO. IS) GO TO 15 
10 CONTINUE 
GO TO 5 

15 WRITE(2fl000) WZ>BL,S.E,IS 
PRINT 1010, WZf8L,StEfIS 

1010 FORMATCIH , FIO . IXt F5 . 3 » IX, E 15 . 8» IX, F 12 . A, IX, 12 ) 

IF (WZ .LT. wn WST0R(1)-WZ 
IF (WZ *GT. (WI ♦CM)) GO TO 5 
L- L+1 
ILINE(L)>tS 
OMEGSTRdl- WZ 
20 IF (WZ.LT.Wn GO TO 5 
K«K ♦! 

WSTOR(K) «WZ 
GO TO 5 

1000 F0RMAT(F10.A,F5.3,E15.8fF12.A,I2) 

25 WRITE(2,1000)WZ,BL,S,E,IS 
RETURN 
END 
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PROGRAM SRTGFCA 


SRTGFCA is the second overlay subordinate to SMART. 

SRTGFCA reads from cards surface temperatures and emissivities/ 
and optical and thermal parameters as required for instrument 
simulation. SRTGFCA reads from temporary storage disk 6 the 
transmittance and radiance calculations performed by TRANSMT, 
and computes instrument response parameters. It is again 
emphasized that any infrared instrument requiring atmospheric 
transmittance and radiance information as input may be sub- 
stituted for SRTGFCA. Storage required is 27342s . 
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□VERLAY(HAIN,2»0I 

PROGRAM SRTGFCA 

DIMENSION ATMTALiC 17, 10),RA0ATM( 17, 10 I , RADSUN < 10# 2 , 10 ) • TAUSL AT ( 10 ) , 
SUNCOMI 10),RADSURFf 100,2,2>#SURFC0M(10, 10>, 

RADT0T( 100 » 2 » 2 I ,SUNW(41>,SUMFLUXC4n,CnNClC10), 

CONCHO), 10(10) 

DIMENSION PLANK(10),BLCKIN(10),VAC8NEW(10),GASBNEW(10) ,BBHNEW(10), 
G AS ANEW (100, 2,2), VAC ANEW! 100, 2, 2), GAS ENEW( 100, 2, 2) , 
BBCNEW( 10),BBH0T(10), BBC0LD1(10),BBHNEWG(10) »BBCNEWG(10) 
,B8HNEWC (10)>BBCNEWC( 10) ,8BHNEWV(10)*BBCNEWV(10), 

BBH0LD( 10),BBC0L0( 10), BBM0LDG(10),BBCnLDG(10) , 

BBHOLDC (10),BBCOLDC(10), VACBOLDdO), GASB01D( 10) , 
GASAOLD(100,2,2),VACAaiD(100,2,2),GASEOlO(100,2,2) , 
RADH0T( 10) ,RADC0LD(10) ,RA0H0TG(10) ,RADC01G(10), 

RADHOTC (10),RADC0LC(10),RADVACB(10),RA0GASB(10), 

RADGASA ( 100,2,2),RA0VACA( 100,2,2),RADGA$E(100,2,2), 
BBHDLDV(l0),BBC0LDV(10),RADH0TV(10),RADCniV(10) 

DIMENSION T0TSUN( 10, 2, 10) , TOTSURF ( 100, 2,2 ) ,RADSUN0( 10, 2, 10), 
T0TRAD( 100,2,2),RADS0LD( 100, 2 , 2 ) , R ADOLO ( 17, 10 ) , 
RADOTOT(100,2,2),TAUOLO( 17, 10), TRANS (17,10 ),T0TATH( 17, 10 
),GBALDIF(10),VBALDIF(10),TAUA(10),SIGBV(10)» 
GCLTAIN(10),ERR0N(10),SIGAV( 100,2,2),0V(100,2,2), 
GCLTAAT(100,2,2),V(100,2,2) 

CaMM0N/AB/TEHP(17) ,PRESC17),THICK(17),0PATH( 17),6ASC0NC(20,17), 
IDENT(20),XMULT( 10), JPR0F(17) 

C0MM0N/AA/NXM,NLAY,NLAYI,N0C,KBI,NSPEC,DELTAW,CMINV,IPOLLUT, 

BROAD, WI,WF 

COMH0M7CA/RESPGt35),FUTERW(35),RESPV(35),RESP(35) 

C0Mn0N/CB/EMlSS(10),TSURF(10) 

DATACSUNW-1941.75, 1960. 20, 2020. 20, 2061. 86, 2105 .26, 2150. 54, 2197*80, 
2247. 19, 2298 .85, 2352 .94, 2409 .64, 2469. 14, 2531 .64, 2597. 40, 
2666. 67, 2739. 73, 2816. 90, 2898 .55, 2985. 07, 3076. 92, 3174 .60, 
3278.69,3389.83,3508.77,3636.36,3773.50,3921.57,4081.63, 
4255.32,4444.44,4651.16,4878.05,5128.21,5405.41,5714.29, 
6060. 61, 6451. 61, 6896. 55, 7407. 41, 8000. 00, 8695. 65) 

DATA(SUNFLUX». 2746-6, .284E-6, .2966-6, • 3Q7E-6, • 323E-6, . 330E-6, 


1 .3436-6,. 3596-6, .3626-6, .3916-6, .4116-6, . 439E-6, 

2 . 4626-6,. 4816-6,. 5016-6,. 5266-6, .554E-6,.580E-6, 

3 .6076-6, .639E-6, .6736-6, .708E-6, .748E-6, ♦ 789E-6, 

4 .8356-6,. 8836-6, .9386-6,. 9956-6, 1.0586-6,1.1286-6, 

5 1.2046-6, I .2866-6, 1.3766-6, 1.4746-6, 1.5816-6, 

6 1.6996-6, 1.8266-6, 1.9646-6, 2. HOE-6, 2. 250E-6, 
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7 2.322E-6) 

DATA C771.1906E-12/.C8/l*439/»C9/6.7675E-8/,Cl0/l09./,NSUN^41/ 

READ 1080»NEMIS»NSURF#N0Cf ZENITH 
1080 F0RHAT(3I5#F10,3) 

READ 1090r ITSURF( J)^ J-l^NSURF) 

1090 F0RHAT(8E10.3) 

READ 1090»f ENISS( J)» J«1>NEMISJ 

READ 1100> (FILTER W(J)#RESPV«J)#RESPG(J)^RESP(JJ»J*1»N0C) 

1100 F0RMAT(Fl0.3>F10.6fF10.6*F10,3) 

READ 1090# IBLCKINCK)»K*1»NIAYI) 

READ 1090# (B8H0T(K)#BBC0LDl(K)#K»l#NLAYn 
READ 1090# EHISBB 

IF (CWl.LT.FILTERWa)).0R.(WF.6T.FILTERW(NDCn> GO TO 500 
JCOUNT-0 
00 5 K-1#NLAY 
DO 5 J-1#NXM 

PAOATMIK# J) -TRANS (K#J)»T0TATM(K#J)»0. 

5 CONTINUE 

00 10 K-l#NLAYI 

RAOHOT(KI*RADCOLO(K|-PAOHOTG(K|.RAOCOLGCK|«RAOHOTC<K>»RAOCOLC(K|. 

1 RAOVACB(K)-RAOGAS0(K)«RADHOTV(KI-RADCOLV(K)*O. 

10 CONTINUE 

DO 15 K-lrNLAY 
DO 15 L-1#NENIS 
00 15 H«1#NSURF 
DO 15 J-1#NXM 
TOTSUN(K#L# J)-0. 

JJ-Jf(K-l)*10 

RADGASA(JJ#L»N)«RADVACA( JJ#L#H)«RADGASE(JJ#L#N)«0. 

TOTRAO( JJ#lf H)*TOTSURF( JJ#L#M)-0. 

15 CONTINUE 
C ♦♦♦ READ DISC 6 ♦♦♦ 

20 READ(6)0HEGA, ( (ATMTAU(Kt J), J-1#NXN)#K-1#NLAY)# ( (RAOATM(K# J)# 

1 J-l#NXm»K-l#NLAYI#DWl 

IF (E0FC6)) 1B5#25 
25 JCOUNT-JCOUNT^l 

C ♦♦♦ INTERPOLATE WAVENUMBER DEPENDENT INPUT VALUES 

CALL FILTERP(OHEGAf TAUVFIL#TAUGFIL#SUNINTP#NOC#NSUN#SUNW#SUNFLUX, 

IRESPIN) 

C ♦♦♦CALCULATE UNATTENUATEO SURFACE EMISSION ♦♦♦ 

CALL PLANKF(OMEGA#SURFCOn#NEHIS#NSURF) 

C ♦♦♦ CALCULATE UNATTENUATEO SOLAR RADIATION REFLECTED BY THE EARTH SURFACE ♦♦♦ 
IF { ZENITH. GE. 90. ) GO TO 30 

CALL SUNNYC OMEGA# NEM IS# SUNCOM#ZENITH#SUNINTP#CHAPZEN> 

GO TO 40 
30 CHAPZEN-0. 

DO 35 L-1#NEHIS 
35 SUNC0M(L)-0. 

C ♦♦♦ CALCULATE TOTAL ATTENUATED RADIATION AT THE TOP OF EACH ATMOSPHERIC 
C LAYER ♦♦♦ 

40 DO 130 K«1#KBI 

DO 130 L-1»NEMIS 
DO 130 M-1#NSURF 
DO 130 J-l»NXM 
JJ-J4(K-1)A10 

TAUSLATC J)-ATMTAU(KBI# J)^^CHAPZEN 

RAOSUN(K#L# J) -SUNCOMCLI ATAUSLATCJ)* ATMTAU(K#J) 

RADSURFt JJ#L#H)«SURFCOM(L»M)^ ATMTAUtK# J J 

RAOTOT(JJ#L»M)-RAOATM(K# J» ♦ R AOSURF C J J# L# M I ♦ R AOSUN (K # L # J » 

130 CONTINUE 

IF (NLAYl.EO.O) GO TO 160 
W3- OMEGA ♦OMEGA ♦OMEGA 

C ♦♦♦ CALCULATE INSTRUMENT DETECTED RADIATION ♦♦♦ 

DO 140 K«l#NLAYt 

PLANKIKI -C7 ♦W3/CEXP(C8 ♦OMEGA /BLCK IN C K 1)-l. ) ♦EMISBB 
VACBNEW(K) -PLANKIKIATAUVEU ♦RESPIN 

GASBNEWtK) -PLANK CK)^ATMTAU(K +KBI#1) ♦TAUGFIL ♦RESPIN 
00 135 L-l»NEMtS 
DO 135 M-1#NSURF 
DO 135 J-1#NXM 
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120 


125 
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1^0 


1^5 


150 


155 


160 


165 


170 


175 


180 


JJ*J ♦^K-l) *10 

GASANEW( JJf t»H)»RAOTOTUJ#L#M) ♦ATMT AUC K+KBI » 1 » ♦TAU6F IL*RESP IN 
VACANEW(JJ9ttN)«RA0T0T(JJ#Lf N) ♦TAUVFIL «RESPTN 
GASENEWC L#N>»(RADT0T(JJ#L»H> *ATNTAU(K^KBT» 1) 4R ADATNCK^KBI* 
1 1) I ♦TAUGFU»RESPIN 

135 CONTINUE 

C ♦** CALCULATE BALANCE AND CALIBRATION PARANETERS *** 

BBHNEW<K}»IC7 ♦W3/<EXP(C8 AOHEGA/BBHOT ( K |)-1 • M ♦TAUVFIL ♦RESPIN 
BBCNEW(K)*<C7 ♦W3/CEXP(CB AONEGA/BBCOLDl ( K H -1 • ) ) ♦TAUVFIL^RESPIN 
BBHNEWGCK)*CC7AW3<<EXPCC8*OHEGA/BBHOTCKn-l. n + ATNTAU<K*KBI#l) + 

1 TAUGFIL+RESPIN 

BBCNEWG(K»«(C7«W3/(EXP(C0«ONEGA/B8COLOl(Kn-l« I )AATNTAU(K4KB I# 1 ) ♦ 

1 TAUGFIL^RESPIN 

BBHNEWC<KI-BBHNEWG(K)*ATHTAUCK+KBI+NLAYI»1)+<RA0ATN(K+KBT+NLAYI,1) 

• ♦TAUGFU*RESPIN)*ATMTAU(K+KBIf 1) 
B8CNEHCIK)*BBCNEWG(K)*ATMTAUCK+KBI+NLAYI>l)f <RADATH<K+KBI+NLAYI»ll 

♦ ♦TAUGFILARESPIN)*ATNTAU(K4^KBI»1) 
BBHNEWV(K}«BBHNEW(K)*ATNTAU(K4-KBI^NLAYI»lK(RA0ATHfK4KBl4NLAYI»l) 

• ♦TAUVFIL*RESPIN> 

BBCNEWV(K)«BBCNEWfK)«ATHTAU(K>fKBI-fNLAYI>l)4(RADATH(K>Kai4>NlAYIf 1) 

♦ ATAUVEIL+RESPIN) 

1^0 CONTINUE 

IF (ONEGA.EOtWn GO TO 1A5 
C ♦♦♦ INTEGRATE INSTRUNENT RADIATION VALUES ♦♦♦ 

CALL XNTE6(0Wl»BBHNEW»BBH0LDf RADHOT»NLAYI ) 

CALL XNTEG(DWl»BBCNEW»BBCOLOf RAOCOLOfNLAYIl 
CALL XNTEG(0Wl>8BMNEWGf BBHOLDG^RAOHOTG^NLAYn 
CALL XNTE6<0Wl^BBCNEWG»BBC0LD6f RAOCOLG^NLAYI) 

CALL XNTEGCDWlf BBHNEWCf BBHOlOCf RAOHOTCfNLAYI) 

CALL XNTEG(DWl»BBCNEWCf BBCOLDC»RAOCOLC#NLAYI) 

CALL XNTEG(OWlf VAC BNEW» VAC BOLDER ADVACB^NL AY n 
CALL XNTEG(DWl#GASBNEW»GASBOLD»RADGASB»NLAYn 
CALL XNTEGCOWl»08HNEWV#BBHOLOV^RAOHOTVf NLAYII 
CALL XNTEG(DWlf B0CNEWV«BBCOLOV#RAOCOLV»NLAYn 
CALL XNTEG2fOWl»GASANEW»GASAOLD» RADGASA^NLAYIfNEWTS^NSURFf NXH) 

CALL XNTEG2 (DWl»VACANEW»VACAOLOf RADVACA»NLAYI»NFNIS»NSURF»NXN) 

CALL XNTEG2C0Wl>GASENEW»GASE0L0»RA0GASEfNLAYl#NEMIS>NSURF»NXH) 

C *♦♦ STORE CURRENT INSTRUMENT RADIATION VALUES FOR INTEGRATION PURPOSES 
1A5 DO 155 K«1^NLAYI 

6BH0LD(K)«BBHNEW(K) 

B6C0L0<K)-BBCNEW«K) 

BBHOLDG(K) -BBHNEWG(K) 

BBCOLDGfK)>BBCNEWG(K| 

BBHOLOC(K)*BBHNEWC(K) 

BBCOLDCCK)*BBCNEWC(K) 

BBHOLDV(K)*BBHNEWV(K) 

BBCOLOV(K) «BBCNEWVCK) 

VACBOLD(K)«VACBNEW(K) 

GAS60LD(K)*GASBNEW(K) 

00 150 L-lfNEHIS 
DO 150 M*1»NSURF 
DO 150 J»lfNXM 
JJ« J+(K-n *10 

GASAOLD(JJ»LfM)»GASANEWUJ»L»N) 

VACAOLO(JJ»L»MI •VACANEWCJJ»L#M) 

GASEOLOfJ J»L» H) «6ASENEW( JJf L*M) 

150 CONTINUE 
155 CONTINUE 

160 IF (OMEGA.EQ.Wn GO TO 165 

C INTEGRATE ATMOSPHERIC RADIATION AND TRANSMITTANCE AND INSTRUMENT 

C TRANSMITTANCE ♦♦♦ 

CALL XNTEG3(DWl»ATMTAU^TAU0L0#TRANSfNLAY»NXH) 

CALL XNTEG3I DWl»RADATM#RAOOLD#TOTATMfNLAY»NXM) 

CALL XNTEG2(0Wl»RA0T0Tf RADOTOTtTOTRADfKBI f NEMIS#NSURF»NXM) 

CALL XNTEG2(DWlfRADSURF»RADS0L0f TQTSURF^KBl » NEMt St NSURF* NXM ) 

CALL XNTEG^(DWl»RAOSUNf RADSUNO#TOTSUN#KBI ^NEMIS^NXM) 

C STORE CURRENT TRANSMISSION AND ATMOSPHERIC RADIATION VALUES FOR 

C INTEGRATION PURPOSES 

165 00 170 K«lfNLAY 
DO 170 J«1>NXM 
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TAUOLD(K»J|*ATHTAU<Kf Jl 
RAOOLD<K»J}«RADAtn(K»JI 

185 170 CONTINUE 

DO 180 K>1»KBI 
DO 175 l«l»NEniS 
DO 175 N-lfNSURF 
DO 175 J»1#MXM 

190 JJ«J^fK-l)A10 

RADOTOTIJJ#l»H)-RAOTOTCJJ#Lf HI 
RADSOLDf JJf L»HI-RAOSURF( 

RADSUNOfK,L#J)«RADSUN(K#L# J) 

175 CONTINUE 

195 180 CONTINUE 

GO TO 20 

C ♦♦♦ CALCULATE INSTRUMENT OUTPUT RESPONSE PARAMETERS 
185 IF (NLAYI.EQ.OI GO TO 200 
DO 195 K.lfNLAYI 

200 GBALDIFCK|«RADHOTG(KI*RAOCOL6CKI 

VBALDIFfKI-RADHOJfKI-RAOCOLOCK) 

TAUA(KI«GBAL0IFCKI/VBALDIF(K) 

SIGBV(KI«RADVACB(KI* tauacki 
GCLTAIN(K|-RADGASB(KI/RADVACB(KI 

205 ERRDN(KI-RAOGASBfK)-TAUA(K)« RADVACBIKI 

DO 190 L*1»MEMIS 
00 190 M»1»NSURF 
Oq 190 J*ltNXM 
JJ«J+ fK-ll* 10 

210 SI6AVUJ#L»M)«RA0VACAUJ#Lf HI* TAUAfKI 

0V( JJfL»M)*RAOGASA(JJ,L,M)- S I GAV ( J L# H ) 
GCLTAATUJ»L»H)«RAOGASACJJ#Lf M)/RADVACAUJ#L^M| 
V(JJ#L#MI-(RAOGASA(JJ#L#MI ♦ S IGA V( J J# L # M 1 1 /?« 

190 CONTINUE 

215 195 CONTINUE 

C *♦♦ OUTPUT 

200 J-0 

00 250 LL«1>16 

IF (IDENTiLLI.EO* 0) GO TO 250 

220 J-J+1 

IF (lOENmil.EQ. 1) IO(J)*2HCO 
IF (IDENTCLLI •EQ. 21 I0UI-3HH20 
IF (IDENTCLLI .EQ. 71 I0(JI*3HN20 
IF HDENTILU .E0.12) IDUI-3HC02 

225 250 CONTINUE 

IF (IP0LLUT.EQ.il IP0LT-2HC0 

IF CIP0LLUT.E0.2I IPOLT*3HH20 

IF (IP0LLUT.E0.3I IP0LT*3H$02 

IF (IP0LLUT.E0.4I rP0LT«3HNN3 

230 IF (IP0LLUT.E0.5) IP0LT»3HCH4 

IF (IP0LLUT.EQ.7I IP0LT-3HN20 

IF (IP0LLUT.EQ.12i IP0LT-3HC02 
C- 2.*CMINV 
OIST-WF-Wl 

235 DO 255 K-1,NLAY 

00 255 J-l#NXM 
TRANS(Kf JI*TRANS(K» JI/OIST 
255 CONTINUE 

PRINT 2000>lP0LTfWIfWFiC#0ELTAW»JC0UNT 

240 2000 F0PMATC1HI#*THE INVESTIGATED POLLUTANT IS AfAS/* THE WAVENUMBER IN 

ITERVAL IS *#F10.3»* TO ♦»F10.3/* THE SUBINTERVAL OF INTEGRATION IS 
2 ♦»Fl0.3f* CM-1*/* THE MINIMUM INTEGRATING INCREMENT IS *fF6*3f* C 
3M-1*/* THE NUMBER OF POINTS OF INTEGRATION IS ♦»I6I 
PRINT 2010 

245 2010 F0RHAT(1H0,*FILTER FUNCTION*/* WAVENUMBER VACUUM RESP GAS CELL 

1 RESP DETECTOR RESP*) 

PRINT 2020» (FILTERW(JIfRESPVUI»RESPGCJ)»RESP(JWJ« 1»N0CI 
2020 FORMATCIH » F 10 . 3» 4X» FIO. 4# 4X# FIO. 4» 6X» FIO. 4 1 
00 300 J«lfNXM 

250 CONCKJI-XMULTm *GASCONC (I POLLUTt 1 1 

300 CONTINUE 
PRINT 2030 
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265 
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280 
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295 


300 


305 


310 


315 


320 


2030 F0R»AT<1H1#*ATM0SPHERIC RARAM6TERS*I 
PRINT 20AO»(IDUIf J«1»NSPEC) 

20A0 FORMATdH >*LAYER TEMP fK> PRESS <ATM» THICK fCM) PATH CATM-CM 
1) A(A3f8X ) } ’ 

DO 310 K-1»NLAY 

BROFAC-CBROAD-l. > ♦GASCONC (I POLLUTE 1 > 

IF (K.LE.KBU PRESf K)«PRES(K) *BR0FAC 
IF (K.GT.KBI) PRES(K)«PRESfK)/BRDFAC 
310 CONTINUE 

DO 330 K-lfNLAY 
J-0 

DO 320 LL«lfl6 

IF (lOENT(LL).EQ. 0) GO TO 320 
J-J+1 

CONCU)-GA$CONC( lOENTlLU^K) 

320 CONTINUE 

PRINT 2050# K» TEMP (K )» PRE S ( K )# THICK (K )» OP ATH( K )>( CONC U )# J -1 # NSPEC I 
2050 FORMATdH , 12# 5X# F7. 2# 3X# F8 . 5# 5X# IPE12 • 5# IX# E12 . 5# 3X# A C ElO, 3# IX » ) 
330 CONTINUE 

PRINT 2060#ZENITH 

2060 F0RMAT(1H0#*THE SUN ZENITH ANGLE IS ♦#0PF7#3#A DEGREES*! 

LOOP*NLAY 

IF (NLAYI.GT.O) LOOP«NLAYI 
DO AOO K*1#L00P 
PRINT 2070#K 

2070 FORMATdHl#*LAYER *#I2» 

IF UPROF(K) .EO.n IPR0F-7HL0RENTZ 
IF UPR0F(K).E0.2) I PROF -5HV0IGT 
PRINT 2075# IPROF 

2075 FORMATdH #*THE SPECTRAL LINE PROFILE IS *#A7! 

PRINT 2080 

2080 FORMATdHO#*ATMOSPHERIC OUTPUT PARAMETERS*) 

PRINT 2090# (CONCIU)# J«1#NXM) 

2090 FORMATdH #*THE CONCENTRATIONS OF THE INVESTIGATED POLLUTANT CONSI 
IDERED ARE*/10X#1P10E12.A) 

PRINT 2100# (TRANS (K#J)#J-1#NXM)# dOTATM(K#J)#J«l#NXM) 

2100 FORMATdH #*TRANS -*# 1P10E12. A# /* TOTATM **#10E12.A) 

DO 3A0 L-1»NEMIS 
PRINT 2110#EMISS(L) 

2110 FORMATdH #*SURFACE EMISSIVITY ■*#0PF6.3) 

PRINT 2120# (TOTSUN(K#L# J)» J-1#NXM) 

2120 FORMATdH #*TOTSUN «*# 1P10E12. A ) 

3A0 CONTINUE 

JJ-1+(K-1)*10 
JK- JJ+9 

DO 350 L-1#NEMIS 
DO 350 M-1#NSURF 
PRINT 2130#EMISS(L)#TSURF(M) 

2130 FORMATdH #*SURFACE EMISSIVITY - *#0PF6.3#/* SURFACE TEMPERATURE - 
1 *#F8.3) 

PRINT 21A0# (TOTSURFf J#L#M)# J-JJ# JKI# (TOTRADf J#L#M|# J-JJ# JK) 

21A0 FORMATdH #*TOTSURF «*# 1P10E12.A# /* TOTRAO -*#10E12.A) 

350 CONTINUE 

IF (NLAYI.EO.O) GO TO AOO 
PRINT 2150 

2150 FORMATdHO#*INSTRUMENT OUTPUT PARAMETERS*) 

PRINT 216O#RAOVACB(K)#RADGASB(K)#TAUA(K)#ERR0N(K) 

2160 FORMATdH #*RA0VACB -*# 1PE12# A# /* RAOGASB ■*#E12.A/* TAUA »*#E 
112 

l.A/* ERRDN ■*#E12.A) 

PRINT 2090»(C0NC1(J)#J«1#NXM) 

PRINT 2170# (TRANS(K4KB1# J)# J«1#NXM) 

2170 FORMATdH #*TRANS •*# 1P10E12. A ) 

DO 360 L-1#NEMIS 
DO 360 M-1#NSURF 
PRINT 2130#EMISS(L)#TSURF(M) 

PRINT ?180#(RADVACA(J#L#M)#J«JJ# JK)#(RADGASA(J#L#M)»J«JJ#JK)i 
1(RADGASE( J#L#M)#J*JJ#JK)# (DV( J#L#M)# J«JJ# JK)# (V(J#L#M)# J«JJ# JK) 
2180 FORMATdH #*RADVACA -*# tP10E12#A/* RADGASA **#10E12#A/* RADGASE «* 
1#10E12.A/* DV •*»10E12«A/* V «*#10E12#A) 

360 CONTINUE 
PRINT 2190 
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325 


330 


335 


3A0 


3<»5 


350 


2190 F0RM*T«IH0»*CAL1BR4TI0N OUTPUT PARAMETERS*) 
OVHOT<RAOHOTC(Kt-TAUA(KI*RAOHOTV(K) 

DVCOLD-RAOCOLC (K ) >TAUA (K ) *RAOCOLV(K ) 

VHOT>(RAOHOTC(K) * TAUA(K)*RAOHOTV(KM /2. 

VC01D-(RADC0LC(K) 4TAUA(KI*RA0C0LVIK)) /2. 

PRINT 2200> TRANS IK*KBl*NLAYI»l)yRADHOTC(K>>RADCOLC(K)>RAOHOTV(K If 
IRAOCOLVIKIf OVHOTf OVCOLOf VHOTf VCOLD 

2200 FORNATdH f*TRANS «*f lPE12.4f /♦ RADHOTC -AfElPtAf/* RADCOLC •*» 
lE12.4f/* RAOHOTV «*fE12.Af/* RAOCOLV <*fE12.9f/* OVHOT >*fE12.Af 
2/A DVCOLD •*fE12.4f/* VHOT •*»E12.4>/* VCOLD •*fE12.4) 

PRINT 2210 

2210 FORMATUH f*BAlANCE OUTPUT PARAMETERS*! 
OVHOTB«RADHOTG«K)-TAUA(K)*RAOHOT«KI 
DVCOLOB«RAOCOL6(K)-TAUA(K)*RAOCOLO(K) 

VH0TB>(RADH0T6(K) ATAUA (K I *R ADHOTI K 1 1 /2 . 

VCni0B*(RA0C0L6(K) ATAUAIK I *R AOCOLO(K 1 1 /2. 

PRINT 2220# RAOHOTlKIfRAOCOlOIKIf RAOHOTGiKIfRAOCOLGIKIfOVHOTBf 
lOVCOLOBfVHOTB# VCOLDB 

2220 FORMATUH #*RAOHOT «*# 1PE12 . 4# /* RAOCOLO «*»E12.4#/* RA0H0T6 •*> 
1E12.4»/* RADCOLG •*»E12.4#/* OVHOTB -*#E12.4#/* OVCOIOR •*»E12.4# 
2/* VHOTB •*#E12.4»/* VCOLDB «*»E12.4) 

400 CONTINUE 
GO TO 550 
500 PRINT 3000 

3000 F0RMATI1H1#*THE WAVENUMBER INTERVAL IS OUTSIDE OF THE FILTER RANGE 
1*1 

550 RETURN 
END 
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SUBROUTINE FILTERP 

FILTER? is a subroutine of SRTGFCA. FILTERP interpolates 
linearly the wavenumber-dependent vacuum cell and gas cell 
filter transmission, detector response, and solar flux values. 
Input filter parameters are read from cards in SRTGFCA while 
solar flux parameters are found in DATA statements in SRTGFCA. 
Storage required is 272s . 


1 SUBROUTINE F ILTERP (W, TAUV> TAUG, SUN,.NOC, NSUN, SUNW,SUNFIUX, RESPIN) 

C ♦♦♦INTERPOLATE WAVENUMBER DEPENDENT VACUUM CEU RESPONSE. GAS CELL RESPONSE. 
C DETECTOR RESPONSE. AND SOLAR FLUX FROH INPUTTED PARAMETERS ♦♦♦ 

DIMENSION SUNW(41I.SUNFLUX(41) 

5 C0MMaN/CA/RESPGt35).FILTERW(35I.RESPV(35).RESP(39t 

DO 10 K-2.N0C 
K2-K 

IF (FUTERW(K).GE.WI GO TO 20 
10 CONTINUE 

10 20 K1-K2-1 

C0H«(W>FILTERWIKin/(FILTERW(K2)-FILTERWIKll I 
RESP1N-RESPIK1I+ COM^ ( RESP (K2 )-RESP (K1 ) I 
TAUV«RESPV(K1) *C0M ♦ < RESPV (K2 )-RESPV( K1 ) ) 

TAUG«RESPG(K1) +C0H ♦ IR ESPG <K2 ) -RESPG t K1 ) ) 

15 IF IW.LT.SUNWIll) GO TO 50 

DO 30 K-2.NSUN 
K2-K 

IF (SUNW(K) .6E.W) GO TO AO 
30 CONTINUE 

20 AO K1.K2-1 

C0M-(W-SUNW(K1)I/(SUNW(K2)>SUNU(K1)) 

SUN>SUNFLUXIK1 ) ♦C0MA(SUNFLUX(K2)-SUNFLUX(K1) ) 

RETURN 
50 SUN-0. 

25 RETURN 

END 
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SUBROUTINE PLANKF 


PLANKF is a subroutine of SRTGFCA. PLANKF calculates 
unattenuated surface emission using Planck's blackbody 
function. The surface temperatures and emissivities are 
read in SRTGFCA. The storage required is 106s • 


1 


5 


10 


SUBROUTINE PLANKF I W»SURFCON>NENISrN$URF) 

C CALCULATE UNATTENUATEO SURFACE EMISSION USING PLANCK FUNCTION 

C nONOCHRONATICALLY 

DIMENSION SURFCON(lOflO) 

C0MM0N/CB/EM1SS(i6).TSURF(I0> 

DATA C7/1.1906E-12/»C8/1.A39/ 

W3«H*HPW 
DO 10 I>I*NEMIS 
DO 10 K>1»NSURP 

SURFC0NII>K>«C7 *W3/ (EXP (C8*W/TSURF (K > I -1 . I 
SUREC0N(1>K)-SURFC0M(I»K) PEMISSItl 
10 CONTINUE 
RETURN 
END 
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SUBROUTINE SUNNY 


SUNNY is a subroutine of SRTGFCA. SUNNY calculates 
attenuated solar radiation reflected by the earth's surface 
using the Chapman function. The sun zenith angle is an input 
to SRTGFCA. The storage required is 608 • 


1 SUBROUTINE SUNNY(W>NEHIS»SUNCOH»rHETA>SUN»CHAPZEM| 

C CALCULATE UMATTEHUATEO SOLAR RADIATION REFLECTED BY THE EARTH SURFACE 

C MONOCHROHATICALLY 

DINENSION SUNCONdOl 

5 CONnONFCB/ENISSaotfTSURFdOl 

COSS -COS dHETA/57. 2957795) 

CHAPZEN- CHAPITHETA) 

DO 10 t-l*NEMIS 

SUNC0N(I)-SUN PCOSS *d.~ENISSd)) 

10 10 CONTINUE 

RETURN 
END 
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FUNCTION CHAP 


CHAP is a function of SRTGFCA called by subroutine- SUNNY. 
CHAP determines slant path transmission by use of the Chapman 
polynomial. The storage required is 53s • 


l FUNCTION CH4P(THETA» 

C ♦♦♦ OETERHINE SLANT PATH TRANSMISSION 

C FOR THETA ,6T. 60. 0EGREES» THE NATURAL LOG OF THE CHAPMAN FUNCTION IS 

C GIVEN 8Y AN 8 TH DEGREE POLYNOMIAL IN THETA 

5 DATA P1/3.1A15<)26535/ 

IF <THETA.6T.60.» GO TO 20 

CHAP-1. 0/C0S<PI *THETA/ieo.l 
RETURN 
20 XX-THETA 

10 P0LY-<m<n-8.39732633E-13*XX +2 . 127A0023E-10I *XX-9. A599A657E-9 ) 

1 *XX-1.37255e04E-6»*XX +3 .30591581E-5 I *XX *2 .30546152E-2 » *XX 

2 -2.AA2119A0IPXX +9. 78869273E+1 )*XX -1. A3635200E+3 
CHAP-EXP(-POLY» 

RETURN 

15 END 
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SUBROUTINES XNTEG^ XNTEG2 , XNTEG3, XNTEG4 


XNTEG, XNTEG2, XNTEG3 , and XNTEG4 are subroutines of 
SRTGFCA. They each integrate monochromatic values using the 
trapezoidal rule, and differ by the dimensions of the 
variables involved in the integration. The storage required 
for XNTEG, XNTEG2, XNTEG3 , and XNTEG4 is 218, 123e, 31 b, and 
77s , respectively. 


1 SUBROUTINE XNTEG< DWl, VALNEW,VALaLO,SUM,K> 

C ♦♦♦ INTEGRATE HONOCHRONATIC VALUES USING TRAPEZOIDAL RULE ♦♦♦ 
DIMENSION VALNEW(10l»VALOLDU0),SUMa0l 
00 50 IK*1#K 

5 SUMIIKI-SUMI IK) f(DWl/2) ♦(VALNEWCIK) •fVALOLDnKM 

50 CONTINUE 
RETURN 
END 


1 SUBROUTINE XNTEG2 CDW1» VALNEW# VALOLDt SUM,KBI »NEMIS,NSURF,NXM ) 

C ♦♦♦ INTEGRATE MONOCHROMATIC VALUES USING TRAPEZOIDAL RULE ♦♦♦ 
DIMENSION VALNEWI 100, 2 ) , VALOLD ( 100, 2, 2 ) , SUM ( 100, 2 ) 

00 50 K«1,KBI 

5 00 50 L*1,NEHIS 

DO 50 M-1,NSURE 
00 50 J-1,NXM 
JJ-J ♦(K-1) *10 

SUMCJJ,L,M) -SUMCJJ,L,M) ♦CDWl/2)* ( VALNEW( JJ,L, M) ♦VALOL0(JJ, 
10 1 L»MI) 

50 CONTINUE 
RETURN 
END 


1 SUBROUTINE XNTEG3 (OWl, VALNE W, VALOLD, SUM,K,NXM) 

C INTEGRATE MONOCHROMATIC VALUES USING TRAPEZOIDAL RULE *** 

DIMENSION VALNEW(17,10),VALQL0(17,10),SUM(17,10) 

DO 50 IK*1,K 

5 DO 50 J«1,NXM 

SUM(IK#J)«SUM(IK,J)+ (DWl/2)* ( VALNEW (IK , J ) -»■ V AL OLD ( IK , J ) ) 

50 CONTINUE 
RETURN 
END 


1 SUBROUTINE X NT E G4 ( DW 1, V ALN E W , VALQLD, SUM , K , N EM I S , NXM ) 

C *** INTEGRATE MONOCHROMATIC VALUES USING TRAPEZOIDAL RULE *** 
DIMENSION VALNEW( 10,2, 10 ) , V ALQLO ( 10, 2 , 1 0 ) , SUM ( 1 0, 2 , 1 0 ) 

00 50 IK»1,K 

5 DO 50 L*1,NEMIS 

00 50 J«1,NXM 

SUM(IK,L, J)*SUM(IK,L,J) * (DWl/2)* ( V AL N E W ( I K , L , J ) + VALOLD(IK,L 
1 , J)) 

50 CONTINUE 

10 RETURN 

END 
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OPERATING INSTRUCTIONS 


Input 

The input from cards is divided into two sections. One 
section is the atmospheric data used by TRANSMT and its 
associated routines. The other section is the instrvunent 
data used by SRTGFCA and its associated routines. 

The physical set-up for the TEMP, PRES, and THICK data 
is arranged such that all atmospheric layers are read in first, 
followed by the instrument cell layers to be considered, 
followed by an equal number of calibration cell layers. An 
option for no instrument and calibration layers is explained 
later. 

In order to avoid unnecessary calculations, one should 
be careful in choosing the values of DELTAW and CMINV. For 
our sample case using carbon monoxide as the primary gas, 
2070-2220 cm“^ as the band width, DELTAW of .01 cm“^ , and 
CMINV of 5 cm~^, a total of 12628 points of integration were 
considered. If DELTAW were chosen smaller and/or CMINV 
were chosen larger this total would be significantly greater 
causing execution time to increase. 

Caution should be exercised concerning the units of the 
line parameters read from the spectral line tape. These units 
must be as specified in the dictionary of FORTRAN variables, 
and the reference temperature of the spectral lines must be 
296. OK or the appropriate change to REFTEMP in the DATA 
statement in TRANSMT must be performed. A sample case input is 
listed in Appendix A. 


^ Options 

Presently, one may choose either Lorentzian or Keilkopf 's 
approximation to the Voigt (ref. 7) profile in the calculation 
of spectral line shape in program TRANSMT. The Doppler line 
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profile will be added in order to more accurately describe 
atmospheric conditions of low pressure and high temperature. 

It should be mentioned here that the Voigt profile requires 
a significant increase in calculation time. 

By setting the value of NLAYI to zero, only atmospheric 
transmittance and radiance calculations are performed. 

The program is currently set up to perform carbon 
monoxide calculations in the 4.6 ym spectral band. Any other 
primary pollutant gas may be considered by setting the value of 
IPOLLUT to the identification nximber of that specie and by 
making necessary adjustments to band widths and filters. The 
gas broadening coefficient, BROAD, must also be changed to 
correspond to the gas being considered. 

Program SRTGFCA simulates a gas filter correlation radio- 
meter. Any infrared instrument using atmospheric transmittance 
and radiance information as input may be substituted for SRTGFCA. 

It should be mentioned here that a routine to plot trans- 
mittance versus wavenumber as a function of altitude may be 
incorporated following the writing of the transmittance values 
on temporary storage disk 6 in TRANSMT. 

The trapezoidal rule is used as the integration approxi- 
mation in routines XNTEG, XNTEG2 , XNTEG3, and XNTEG4 . Should 
the user wish to incorporate a different integration process, 
these routines may be easily replaced. Similarly, that sub- 
routine that determines the wavenumber mesh for integration 
may be replaced by an algorithm more suitable to the user . 

Output 

The program output is in two sections. The first section 
lists the spectral lines used in the calculations, and the 
input data such as atmospheric profile and filter parameters. 

The second section is divided into two parts. Atmospheric 
output parameters, i.e., band integrated transmission and 
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total upwelling radiance as well as the band integrated com- 
ponents of the total radiance (solar, surface, and atmospheric) , 
as a function of surface temperatures and emissivities are 
listed for the top of each atmospheric layer. The second part 
lists the corresponding instrxment output responses in addition 
to calibration and balance output parameters as required by 
input instructions. If one sets NLAYI to zero, the instrument 
calibration will not be performed. 

A sample case output is listed in Appendix B. 

SAMPLE PROBLEM 


The input listed in Appendix A and output in Appendix B 
correspond to a problem of determining the GFCR response of the 
instrument flying at the top of each of the first four modeled 
atmospheric layers. The model chosen is a 45® North Latitude 
July atmosphere with a corresponding water vapor profile 
(ref. 10) . The sun zenith angle is 45®. The pollutant is 
carbon monoxide calculated at 10 different concentrations in 
each layer, with interferent constituent concentrations of 
water vapor and carbon dioxide held constant in each layer. 

The spectral band is 4.6 ym, and the filter parameters simulate 
those of a NASA/LaRC version of a GFCR. 

COMPARISONS AND CONCLUSIONS 


Obviously, any computer program is only as accurate as 
the theoretical model on which it is based and the accuracy 
of the niimerical algorithms coded into program instructions. 
As outlined previously in the Problem Definition and Method 
of Solution section, the accuracy of this model is dictated 
by the assumptions made in our application to a specific 
problem,. 

In order to verify the results of the SMART calculated 
absorption coefficients, several comparisons were made with 
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existing experimental and theoretical absorption results in 
the CO fundamental band. The SMART program's algorithms for 
computing transmittance were checked against hand calculations 
for several wavenumbers and were shown to be correct. The 
high resolution absorption measurements, .05 cm~^ , reported by 
Chaney and Drayson (refs. 11 and 12) for the R20 line of the 
4.6 ym CO band was simulated by the SMART program, degraded, 
and compared. The absorption was computed at .01 cm~^ intervals 
using the AFCRL line parameter data (ref. 5) , then convoluted 
with a triangular slit transmission function whose width was 
.12 cm~^ at half-maximum. The discrepancy found between the 
higher values of line strengths by Drayson (ref. 12) which were 
calculated from the measurements of Chaney (ref. 11) , and the 
lower line strengths from the AFCRL line parameter tape, i.e., 

1.11 and 0.80 (atm“^cm“^ ) g^p, respectively, for the R20 line, 
will be investigated further in the future. A preliminary 
comparison has shown no discrepancy between Drayson 's line 
strengths near the CO fundamental band center, e.g. , the RO 
and PI line. As the distance of a line from the band center 
increases, the discrepancy increases between the higher values of 
Drayson 's line strengths and the lower AFCRL data line strengths 
(refs. 13, 14, 15, and 16) . The SMART calculated spectrum is 
in excellent agreement with the experimental spectrvim near the 
line center, and has only slightly higher values of absorption 
in the line wings. This deviation in the wings is probably due 
to the difference between the actual line profile (ref. 17) and 
the Lorentz line profile employed in our calculations. The 
absorption in the line wings calculated by SMART is in agreement 
with Burch and Gryvnak (ref. 18) , who have shown that the 
extreme wings of the spectral lines are sub-Lorentzian. Despite 
the line wing shape differences, a comparison between the 
equivalent line width calculation using the Ladenberg and Reiche 
function (ref. 19) as described by Kondrat'Yev (ref. 20) , and 
the SMART integrated absorption reveals less than 0.5 percent 
discrepancy in the total absorption by the R20 line under Chaney's 
test conditions. 
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SMART calculations of the entire CO fundamental band were 
performed and compared to the integrated absorption reported 
by the experimental work of Burch and Gryvnak (ref. 21) . The 
SMART integrated absorption differed from the value reported by 
Burch by only 6.0 percent. This difference is probably due to 
a combination of an uncertainty in the reported experimental 
conditions and Burch's reported ±5.0 percent instrumentation 
measurement error (ref. 21) . 

A line-by-line radiative transfer program that simulates 
a gas filter correlation radiometer has been developed and 
described. The overlaid structure allows for the substitution 
of the instrument simulation program with one that describes any 
infrared sensor. Any of the specific task subroutines may 
easily be substituted by algorithms that are more suitable and 
convenient to the user. 

The program is currently constructed to perform carbon 
monoxide calculations in the fundamental 4.6 ym spectral band, but 
by making necessary adjustments, other gas species transmittance 
calculations may be made. Future additions to the program will 
include the Doppler profile option (Lorentz and Voigt are 
presently incorporated) , a sub-Lorentzian wing absorption al- 
gorithm, the water vapor and nitrogen continuum absorption 
algorithms, and a routine to plot transmittance as a function of 
wavenumber . 

The storage required for the program is dependent upon the 
number of layers (atmospheric, instriiment cell, and calibration 
cell) desired for consideration. For our sample case, 17 layers 
required 130300 octal words of storage and 838 seconds of 
execution time. Although this may appear to be a large amount 
of time and storage, the present program is equivalent to 32 
separate computer runs of an earlier version of a radiative 
transfer program used at LaRC. By simultaneously performing 
several functions, such as balancing and calibrating, program 
efficiency increased, and a cost savings factor of five was 
realized. 
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A sample problem was processed. The resulting computer 
listing of input and output is shown. 
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APPENDIX A 


Input Parameters for a Sample Test Case Listed as Card Images 


17 4 31 10 3 12 12 

1 1 I 1 1 1 I 1 1 1 I I I 1 1 I I 
2C7C. 222a. .01 5. 1 1.00 

0. .1 .2 .3 .a .5 ,7 

• 0 1.0 

293.5 .9831 3C4fl0. 

292.1 .9497 33400. 

290.7 .9163 30480. 

289.3 .883“^ 30400. 

288.0 .8532 30480. 

286.6 .8228 30400. 

205.1 .7927 30480. 

283. 4 .7649 30400. 

26 I .6 .7375 30400. 

31 I .5 .92105 1 . 

31 1 .5 .92105 1 . 

31 1 .5 .9210® 1 • 

31 1 .5 .92105 1 . 

290.5 .92105 2. 

290.5 .92105 2. 

296.5 .92105 2. 

298.5 .92105 2* 

1.0E-C6 1.S44E-C2 

l.Ce-06 1.676E-C2 

3.2E-04 

1 .0£-06 I .498E-02 

3.2E-04 

1 .CE-C6 1 .331E-02 
I .C=:-C6 I .21 3E-02 

3.2E-04 

1.0E-C6 1.092E-02 

3.2E-04 

I .OE-06 9.6755-03 

3.2E-04 

l.OE-06 8.606E-C3 

3.2E-04 

1 .0E-C6 7.520E-C3 

3.2E-04 

.35 

.35 
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• 35 


• 35 

• 08143 

• 08143 


• 08143 

• 08143 


2 2 31 

45 

• 

300. 296. 



•98 .88 



2070.0 

0.0 

0,0 

2075.0 

.0050 

• 0C5C 

2080.0 

• 025C 

• 0250 

2085.0 

• 1 060 

• 1 060 

2090.0 

. lOOC 

. 1 000 

2095.0 

.1550 

• 1 550 

21 OO.O 

.2300 

.2300 

2105.0 

.3500 

• 3500 

2110.0 

.4850 

.4050 

2115.^ 

.6250 

.6250 

2120.^ 

• 7750 

• 7750 

2125.0 

• 8700 

.8700 

2130.0 

.9300 

.9300 

2135.0 

.9700 

• 9700 

2 1 4 0 . n 

.9950 

• 9950 

2145.0 

1 .000 

1 .C?0 

2150. o 

.^950 

• ^5950 

2155.0 

• 9 6 5 C 

• 9650 

2 1 6C .0 

.9050 

• 9C5C 

2 1 65.0 

.8100 

• 8130 

2l7r .0 

.6800 

• 6800 

2175.0 

.5350 

• 535n 

2180.0 

.4 1 OC 

.4 1 OC 

2165.0 

.2850 

.2850 

2 1 50.0 

• 1 9C0 

. 1 900 

2195.C 

. laOO 

.1400 

2200.0 

• C950 

.0950 

2205. C 

.C65 ■'1 

.0650 

2210. 0 

.0350 

• 0330 

2215.0 

.0 150 

• 0150 

2220.0 

• 0050 

• 0C5D 

31 1.3 311 .8 

31 1 .8 

308.8 2e4•^ 

' 308.8 

1 • 




I 

I 

I 

1 

1 

I 

1 

1 

1 

1 

1 

I 

1 

1 

1 

I 

1 

I 

I 

1 

I 

I 

I 

1 

I 

1 

1 

1 


t • 

1 • 

8 

7 30a«8 284.7 3C8.8 


284 , 
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APPENDIX B 


Output Listing for a Sample Test Case 


THE INVESTIGATED POLLUTANT IS CO 

THE WAVENUMBER INTERVAL IS 2070.000 TO 2220.000 
THE SUBINTERVAL OF INTEGRATION IS 10.000 CH-l 

THE NINIMUM INTEGRATING INCREMENT IS .010 CM-l 
THE NUMBER OF POINTS OF INTEGRATION IS 12628 


FILTER FUNCTION 


WAVPNUHBER 

VACUUM RESP 

GAS CELL RESP 

DETECTOR RE 

2070.000 

0.0000 

0.0000 

1.0000 

2075.000 

.0050 

• 0050 

1.0000 

2080.000 

.0250 

.0250 

1.0000 

2085.000 

.1060 

• 1060 

1.0000 

2090.000 

.1000 

.1000 

1.0000 

2095.000 

.1550 

.1550 

1.0000 

2100.000 

.2300 

• 2300 

1.0000 

2105.000 

.3500 

.3500 

1.0000 

2110.000 

.4850 

.4850 

1.0000 

2115.000 

.6250 

.6250 

1.0000 

2120.000 

.7750 

.7750 

1.0000 

2125.000 

.8700 

.8700 

1.0000 

2130.000 

.9300 

.9300 

1.0000 

2135.000 

.9700 

.9700 

1.0000 

21^0.000 

.9950 

.9950 

1.0000 

2145.000 

1.0000 

1.0000 

1.0000 

2150.000 

.9950 

• 9950 

1.0000 

2155.000 

.9650 

.9650 

1.0000 

2160.000 

.9050 

.9050 

1.0000 

2165.000 

.8100 

.8100 

1.0000 

2170.000 

.6800 

.6600 

1.0000 

2175.000 

.5350 

.5350 

1.0000 

2180.000 

.4100 

.4100 

1.0000 

2185.000 

.2850 

.2850 

1.0000 

2190.000 

.1900 

.1900 

1.0000 

2195.000 

.1400 

.1400 

1.0000 

2200.000 

.0950 

.0950 

1.0000 

2205.000 

.0650 

.0650 

1.0000 

2210.000 

.0350 

.0350 

1.0000 

2215.000 

.0150 

.0150 

1.0000 

2220.000 

.0050 

.0050 

1.0000 
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ATMOSPHERIC PARAMETERS 


LAYER 

TEMP (KI 

PRESS (ATM) 

THICK (CM) 1 

PATH (ATH-CM) 

CO 

H20 

C02 

1 

293*50 

.98310 

3.0A800E4>04 

2.996A9E40A 

l.OOOE-06 

l.BAAE-02 

3.200E-0A 

2 

292.10 

•9A970 

3.0A800E40A 

2.89A69E40A 

1.000E-06 

1.676E-02 

3.200E-0A 

3 

290.70 

.91630 

3.0A800E40A 

2.79208E4OA 

1.000E-06 

1.A98E-02 

3.200E-0A 


209.30 

.08370 

3.0A800E40A 

2.69352E+0A 

l.OOOE-06 

1.331E-02 

3.200E-0A 

5 

288.00 

.85320 

3.0A800E40A 

2.60055E40A 

l.OOOE-06 

1.213E-02 

3.200E-0A 

6 

286.60 

•82280 

3.0A800E40A 

2.50789E+04 

l.OOOE-06 

1.092E-02 

3.200E-0A 

7 

285.10 

.79270 

3.0A800EAOA 

2.41615E40A 

1.000E-06 

9.675E-03 

3.200E-0A 

8 

283. AO 

.76A90 

3.0A800E40A 

2.331A2E40A 

l.OOOE-06 

8.608E-03 

3.200E-0A 

9 

281.60 

.73750 

3.0A800E40A 

2.2A790E40A 

1.000E-06 

7.528E-03 

3.200E-0A 

10 

311.50 

.92105 

1.00000E400 

9.21050E-01 

3.500E-01 

0. 

0. 

11 

311.50 

.92105 

1.00000E400 

9.21050E-01 

3.500E-01 

0. 

0. 

12 

311.50 

.92105 

1.00000E400 

9.21050E-01 

3.500E-01 

0. 

0. 

13 

311.50 

•92105 

1.00000E400 

9.21050E-01 

3.500E-01 

0. 

0. 

1^ 

298.50 

.92105 

2.00000E400 

1.8A210E400 

8.1A3E-02 

0. 

0. 

15 

298.50 

.92105 

2.00000E400 

1.8A210E400 

8.1A3E-0Z 

0. 

0. 

16 

298.50 

.92105 

2.0000DE400 

1.8A210E400 

8.1A3E-02 

0. 

0. 

17 

298.50 

.92105 

2.00000E400 

1.0A21OE4OO 

0.1A3E-O2 

0. 

0. 


THE SUN ZENITH ANGLE IS A5.000 DEGREES 


57 



LAYER I 

THE SPECTRAL LINE PROFILE IS lORENTZ 


ATHOSPHERIC OUTPUT PARAMETERS 

THE CONCENTRATIONS OF THE INVESTIGATED POLLUTANT CONSIDERED ARE 



0. 

l.OOOOE-07 

2.0000E-07 

3.0000E-07 

<».0000E-07 

5.0000E-07 

6.0000E-07 

7.0000E-07 

0.OOOOE-O7 

l»0000E-06 

TRANS 

• 9.0869E-01 

9,0950E-01 

9.0053E-01 

8.9676E-01 

8.9318E-01 

8^8976E-01 

0.865OE-OI 

8,0337E-O1 

e^803eE->OI 

8^7i»76E-01 

TOTATM 

SURFACE 

• ^.7195E-06 

EMISSIVITY • 

^•9102E-O6 
• 980 

5.1066E-06 

5.2856E-06 

5.A559E*06 

5-6183E-06 

5.7733E-06 

5.9217E-06 

6^0638E-06 

6.331<tE-06 

TOTSUN 

SURFACE 

■ ^.8823E-07 

EMISSIVITY » 

4.6696E-07 
• 880 

^. 538<»E-07 

<t.<i<fO^E-07 

^.3596E-07 

^•2897E-07 

^•2273E-07 

^.1705E-07 

^•neoE-07 

<».0Z31E-07 

TOTSUN 

SURFACE 

SURFACE 

- 2.929AE-06 

EMISSIVITY • 
TEMPERATURE • 

2.8017E-06 
• 980 
SOO^OOO 

2^7230E-06 

2.66^2E-06 

2^6158E-06 

2^5730F-O6 

2^536^E-06 

2.5023E-06 

Z.^708E-06 

2-^139E-06 

TOTSURF 

■ 5.3556E-05 

5.3311E-05 

5^3079E-05 

5^2859E-05 

5.26<r9E-05 

5^24*i9E-05 

5.2259E-05 

5.2076E-05 

5.1901E-05 

5.1572E-05 

TOTRAO 

surface 

SURFACE 

» 5.8763E-05 

EMISSIVITY • 
TEMPERATURE ■ 

5^8696E-05 
• 980 
296^000 

5.86^0E-05 

5.B588E-05 

5.8S<»lE-05 

5^0f97E-O5 

5^8^55E-05 

5^0^15E-O5 

5^8377E-05 

5^8306E-05 

TOTSURF 

• ^».66l7E-05 

^.6^0^E-05 

^.6202E-05 

<t.6010E-05 

^.58Z0E-O5 

^•565^E-05 

^•5^t88E-05 

^•5329E-05 

A*5177E-05 

^•4890E-05 

TOTRAD 

SURFACE 

SURFACE 

« 5*102^E-O5 

EMISSIVITY • 
TEMPERATURE - 

5.1769E-05 
• 880 
300^000 

5^1763E-05 

5^17A0E-05 

5*1720E-05 

5.1701E-05 

5.168AE-05 

5.1660E-O5 

5.1653E-05 

5^1624E-05 

TOTSURF 

It ^.B091E-05 

^•7871E-05 

^•7663E-05 

^•7^65E-05 

^•7277E-05 

^.7097E-05 

^.6926E-05 

^.6762E-05 

^•6605E-05 

A.6309E-05 

TOTRAD 

SURFACE 

SURFACE 

- 5.57A0E-05 

EMISSIVITY • 
TEMPERATURE • 

5.5591E-05 
• 680 
?96^000 

5.5493E-05 

5.5A15E-05 

5»5399E-05 

5.5290E-05 

5.5236E-05 

5.5186E-05 

5.51A0E-05 

5.5055E-05 

TOTSURF 

• ^•l860E-05 

A.1669E-05 

^•Ue8E-05 

A.1315E-05 

^•U52E-05 

<»^0995E-05 

A.08<»6E-05 

^•070^E-05 

^.0567E-05 

<|.0310E-05 

TOTRAD 

• A.9509E-05 

^•9389E-05 

^•9317E-05 

4.9265E-05 

^•9223E-05 

««9100F-O5 

^•9156E-05 

A.9128E-05 

^.91026-05 

^•9055E-05 


INSTRUMENT OUTPUT PARAMETERS 
RADVACB • ^.0B2AE-09 

RADGASB - 3.3833E-05 

TAUA ■ 8.?069E-O1 

ERRON . 2.6392E-09 

THE CONCENTRATIONS OF THE INVESTIGATED POLLUTANT CONSIDERED ARE 




0. 

l.OOOOE-07 

Z.OOOOE-07 

3.0000E-07 

A.OOOOE-07 

5,0000E-07 

6.0000E-07 

7.0000E-07 

8.0000E-07 

1.00006-06 

TRANS 

■ 

0.^2Z6E-Ol 

e.42Z6E-01 

0.<»ZZ6E-Ol 

0.AZ26E-O1 

8.92Z6E-01 

8.A226E-01 

8.^226E-01 

8.4226E-01 

e.A226e-01 

8.42266-01 

SURFACE 

EMISSIVITY • 

• 960 









SURFACE 

TEMPERATURE ■ 

300.000 









RAOVACA 

m 

2.6875E-05 

Z.60AOE-O5 

Z.6810E-05 

2.6704E-O5 

Z.6759E-05 

2.6736E-05 

2.6714E-05 

2.6693E-05 

2.6674E-05 

2.66376-05 

RAD6ASA 

s 

Z^2270E-05 

Z.Z263E-09 

Z.ZZ56E-05 

2.2Z49E-05 

2.22^36-05 

2.2236E-05 

2.2230E-05 

2.22Z3E-05 

2.2217E-05 

2.22056-05 

RAOGASE 

■ 

2.9195E-09 

2.9187E-05 

Z.9180E-05 

2,917AE-05 

2.9167E-05 

2.9T61E-05 

2.915^6-05 

2.91<»0E-O5 

2. 91426-05 

2.9129E-05 

DV 

m 

-1.Z146E-09 

2.0897E-08 

3.0^68E-O8 

5.3762E-08 

6.75^8E-0a 

8.0173E-08 

9.1835E-08 

l,0267E-07 

1.1278E-07 

1.31116-07 

V 

m 

Z^ZZ71E«09 

Z.Z253E-05 

2.2237E-05 

2.2Z2ZE-05 

2.2209E-05 

2.219ftP-05 

2.218<»E-05 

2.2172E-05 

Z.2161E-05 

2.21396-05 

SURFACE 

EMISSIVITY - 

• 900 









SURFACE 

TEMPERATURE ■ 

Z96«000 









RAOVACA 

■ 

Z.3686E-05 

Z.3666E-09 

Z.365^E-05 

2.36<»2E-05 

2.3632E-05 

2.3622E-05 

2.3613E-05 

2.3605E-05 

2.3597E-05 

2.35826-05 

RADGASA 

■ 

1.9630E-05 

1.9625E-05 

1.9622E-05 

1.961BE-05 

1.9615E-05 

1.9ftllE-05 

1.9608E-05 

1.9605E-05 

1.9602E-05 

1.95956-05 

RAOGASE 

a 

2.6554E-05 

Z.6550E-05 

2.65<»6E-05 

2.65A3E-05 

2.6539E-05 

2.653ftF-05 

2.6532E-05 

2.6529E-05 

2.6526E-05 

2.65Z0E-05 



DV 

■ 

1.0087E-09 

142903E-08 

2.0183E-0B 

2.62A7E-0B 

V 

« 

1.9629E-05 

l*9619E-05 

1.9612E-09 

1.9605E-09 

SURFACE 

EMISSIVITY • 

«6B0 



SURFACE 

TEMPERATURE • 

300.000 



RADVACA 

• 

2.5501E-05 

2.5A22E-05 

2.9371E-09 

2.5331E-05 

RAOGASA 

a 

2.1127E-05 

2.1110E-05 

2.109^E-05 

2.1078E-05 

RAOGASE 

• 

2.8092E-09 

2aB03AE-05 

2.8018E-09 

2.8003E-05 

OV 

■ 

-A.B932E-09 

A.27A7E-08 

6.8991E-08 

8.6792E-08 

V 

■ 

2.1130E-09 

2a089E-05 

2.1059E-05 

2.1039E-09 

SURFACE 

EHISSIVITY • 

.880 



SURFACE 

TEMPE»*TURE • 

296.000 



RAOVACA 

« 

2.2637E-05 

2.257AE-05 

2.2536E-09 

2.2510E-05 

RAOGASA 

• 

1.B756E-05 

1.87A2E-05 

1.8728E-05 

1.8716E-09 

RAOGASE 

• 

2.9681E-05 

2 .96666-09 

2.9693E-05 

2.56A0E-09 

OV 

■ 

-2.B968E-09 

3.9210E-08 

5.2932E-08 

6.208AE-08 

V 

• 

1.8758E-05 

1.872AE-09 

1.8702E-05 

1.8605E-O5 


CALIBRATION OUTPUT PARAMETERS 
Cji TRANS • B.9319E-OI 

KD RAOHOTC • 3.0296E-05 

RAOCOLC • 1.3627E-05 

RAOHOTV • 3.5813E-05 

RADCOLV • 1.7191E-09 

OVHOT • 6.1B33E-07 

OVCniO « -5»B639E-07 
VHOT ■ 2.99B7E-05 

VCOLD • I.3920E-05 

BALANCE OUTPUT PARAMETERS 
RAOHOT ■ 3.7086E-05 

RAOCOLO • 1.3933E-09 

RAOHOTG ■ 3.0736E-05 

RAOCOLG - 1.3206E-05 

OVHOTB • 2.826BE-09 

OVCQLDB • 2.826BE-09 

VHOTB ■ 3.073AE-05 

VCOL08 • 1.3205E-05 


3.1<i05E-0a 3.5952E-08 

I.9599E-05 1.9593E-09 


2.5296E-05 2.5266E-09 
2.1063E-05 2.10A9E-05 
2.7988E-05 2.7973E-05 
l.OOAlE-07 1.11A6E-07 
2.1013E-09 2.0993E-05 


2*2^S8E-05 2.2A69E-09 
1.870AE-05 1.8692E-05 
2.5628E-05 2.5617E-05 
6.7951E-08 7.1757E-OB 
1.8670E-05 1.8656E-05 


A.00A2E-08 ^.3766E-0B 

1.9588E-05 1.9583Er05 


2*5238E-05 2.5212E-05 
2.1035E-05 2.1022E-05 
2.7960E-05 2.79^6E-05 
1.2079E-07 1.288AE-07 
2.0975E-05 2.0958E-05 


2.2A53E-05 2.2A39E-05 
ltB681E-05 1.6671E-05 
2.5606E-05 2.5595E-05 
7.A260E-08 7,99ABE-0B 
1.B6AAE-05 1.0633E-O5 


A.7184E-Oa 5.3273E-08 

1.9578E-05 1.9569E-05 


2.51B8E-05 2.5UAE-05 
2.1009E-05 2.0989E-05 
2.793AE-05 2.7909E-05 
1.3392E-07 l,A786E-07 
2*09AlE-05 2,0911E-05 


2.2A29E-09 2.2A0IE-05 
1.B661E-05 1*66A1E*09 
2.5565E-05 2.5566E-05 
7.7O10E-O8 7. 79696-08 
1.B622E-05 1.8602E-09 


LAYER 2 

THE spectral LINE PROFILE IS LORENTZ 


ATMOSPHERIC OUTPUT PARAMETERS 

THE CONCENTRATIONS OF THE INVESTIGATED POLLUTANT CONSIDERED ARE 



0. 

l.OOOOE-07 

2.0000E-07 

3.0000E-07 

9.0000E-07 

5.0000E-07 

6.0000E-07 

7.0000E-07 

8.0000E-07 

1.00006-06 

TRAMS 

■ S. 65^06-01 

8.5771E-01 

8.5079E-01 

8.9450E-01 

0.3875E-O1 

8.33A6E-01 

8.2856E-01 

0.2AOOE-OI 

0.19726-01 

8.1190E-01 

TOTATM 

SURFACE 

■ 6.7130E-06 

EMISSIVITY ■ 

7.0666E-06 
• 980 

7.3850E-06 

7.6730E-O6 

7.9377E-06 

8.180AE-06 

8.A051E-06 

8.61A2E-06 

8.8099E-06 

9.16786-06 

TOTSUN 

SURFACE 

• A.77^3E-07 

EMISSIVITY • 

A.5535E-07 
• 880 

A.A200E-07 

A.3209E-07 

9.2307E-O7 

%. 16776^07 

A.10A3E-07 

A.0A66E-07 

3.9939E-07 

3.09706-07 

TOTSUN 

SURFACE 

SURFACE 

• 2.86A6E-05 

EMISSIVITY • 
TEMPERATURE ■ 

2.7321E-06 

.980 

300.000 

2.6520E*06 

2.5923E-06 

2.9A32E-06 

2.5006E-06 

2.A626E-06 

2.^280E-06 

2.3960E-06 

2.3302E-O6 

TOTSURF 

• 5.0853E-05 

3.0A06E-09 

5.0003E-05 

4.9637E-05 

A,9302E-05 

A.8993E-09 

A.8708E-05 

A.8AA1E-05 

9.8192E-05 

9.7736E-05 

TOTRAD 

SURFACE 

SURFACE 

• 5.80A3E-05 

EMISSIVITY • 
TEMPERATURE • 

5.7928E-05 
• 980 
296.000 

9.7830E-09 

9.77AZE-09 

5.7663E-05 

5.7990E-05 

5.7523E-05 

5.7960E-05 

9.7901E-05 

5.72936-05 

TOTSURF 

• A.A262E-05 

A.3B73E-05 

A.3522E-09 

A.3203E-09 

A,2912E-05 

A.26A3E-05 

A.2395E-05 

A.2163E-05 

9.1996E-05 

9.15996-05 

TOTRAD 

SURFACE 

SURFACE 

• 5.1A52E-05 

EMISSIVITY • 
TEMPERATURE • 

3.139AE-09 
• 880 
300.000 

9.13A9E*09 

5.13096-09 

5,1273E-09 

5.12A1E-09 

5.1210E-05 

5.1182E-09 

5.11566-05 

5.11076-05 

TOTSURF 

« A.56GAE-05 

A.5262E-05 

A.A900E-09 

A.A572E-09 

9.A271E-05 

A.399AE-05 

9.3737E-09 

9.39906^05 

9.32756-05 

9.2865E-05 

TOTRAD 

SURFACE 

SURFACE 

• 5.52AIE-05 

EMISSIVITY ■ 
TEMPERATURE • 

5.5061E-05 

.880 

296.000 

5.A937E*05 

5.AB38E-05 

5.A752E-09 

5.A675E-05 

5.9605E-05 

5.9591E-09 

5.99806-05 

5.937IE-05 

TOTSURF 

• 3.97^5E-05 

3.9396E-05 

3.9081E-09 

3.8795E-05 

3.8533E-05 

3.B292E-09 

3.8069E-09 

3.7861E-09 

3.7666E-05 

3.73096-05 

TOTRAD 

- 4.9323E-05 

9.919AE-09 

A.9118E-09 

9.9061E-05 

A.901AE-05 

A.8973F-05 

9.0936E-O5 

9. 09036-05 

9.8872E-05 

9.00156-05 


O INSTRUMENT OUTPUT PARAMETERS 

RADVACB « A.082AE-05 

RAD6ASB - 3.3B33E-09 

TAUA - 8.2a69E-01 

ERRDN • 2.6392E-09 

THE CONCENTRATIONS OF THE INVESTIGATED POLLUTANT CONSIDERED ARE 


TRANS 

■ 

0. 

8.9226E-01 

l.OOOOE-07 

8.9226E-01 

2.0000E-07 

8.9226E-01 

3.0000E-07 

8.9226E-01 

9.0000E-07 

8.9226E-01 

5.0000E-07 

B.9226E-01 

6.00006-07 

8.9226E-01 

7.0000E-07 

8.92266-01 

8.00006-07 

8.9226E-01 

1.00006-06 

8.9226E-01 

SURFACE 

SURFACE 

RADVACA 

EMISSIVITY • 
TEMPERATURE • 
• 2.6562E-05 

• 900 

300.000 

2.65016-05 

2.69506-05 

2.6909E-05 

2.63636-05 

2.6325E-05 

2.6290E-05 

2.6257E-05 

2.6227E-05 

2.6171E-05 

radgasa 

■ 

2.2007E-05 

2.19956-05 

2.1903E-O5 

2.1971E-05 

2.1960E-05 

2. 19996-05 

2.1930E-O5 

2.1927E-05 

2.19166-05 

2.18956-05 

RADGASE 

■ 

2.8931E-05 

2.89196-05 

2.09076-05 

2.08966-05 

2.08896-05 

2.88736-05 

2.0862E-O5 

2. 08516-05 

2.88916-05 

2.88206-05 

DV 

• 

-9.8006E-09 

3.35296-OB 

6.92B3E-08 

9.09686-08 

1.13306-07 

1.3397E-07 

1.5199E-07 

1.6758E-07 

1.02156-07 

2. 07966-07 

V 

« 

2.2009E-05 

2.1978E-05 

2.1951E-05 

2. 19266-05 

2.1903E-05 

2.1802E-O5 

2.1062E-O5 

2.18936-05 

2.18256-05 

2.1792E-05 

SURFACE 

SURFACE 

RADVACA 

EMISSIVITY • 
TEMPERATURE ■ 
- 2.3529E-05 

.980 

296.000 

2.3999E-05 

2.3970E-05 

2.39996-05 

2.39306-05 

2.3913E-05 

2.3398E-05 

2.33836-05 

2.33696-05 

2.3399E-05 

RADGASA 

■ 

1.99996-05 

1.99876-05 

1.99816-05 

1.99756-05 

1.9970E-05 

1.9969E-05 

1.9959E-05 

1.9953E-05 

1.9998E-05 

1.99386-05 

RADGASE 

■ 

2.69186-05 

2.69126-05 

2.6906E-05 

2.6900E-05 

2.6399E-05 

2.6389E-05 

2.6383E-05 

2.63786-05 

2. 63736-05 

2.6363E-05 



DV 

« 

-9.1766E-10 

1.B3A0E-08 

3.2138E-09 

^.3^63E-08 

V 

a 

1.9^9AE-05 

1.9478E-05 

1.9969E-09 

1.9S5AE-09 

SURFACE 

EMISSIVITY ■ 

• 880 



SURFACE 

TE»1PER*TURE • 

300.000 



RADVACA 

8 

2.5283E-05 

2.51B8E-05 

2.9123E-05 

2. 50716-05 

RADGASA 

a 

2.09A5E-05 

2.092AE-05 

2.090AE-09 

2.08856-05 

RAOGASE 

a 

2.7869E-05 

2,78A8E-09 

2.78Z9E-05 

2.701OE-O5 

OV 

a 

-6.8533E-09 

9.1010E-0B 

B.A677E-08 

1.0872E-07 

V 

a 

Z.09A8E-05 

2.0898E-09 

2.0B62E-05 

2.0831E-05 

SURFACE 

EMISSIVITY ■ 

• 880 



SURFACE 

TEMPERATURE • 

296.000 



RADVACA 

a 

2.2595E-05 

2.2AB7E-09 

2.2AA7E-09 

2.2^10E-O5 

RADGASA 

a 

1.B688E-05 

1.B672E-09 

1.8658E-05 

l.B64^E-09 

RADGASE 

a 

2.5612E-05 

2.9997E-05 

2.99B2E-05 

2.5568E-09 

OV 

a 

-3.007AE-09 

3.7371E-08 

9.5812E-08 

6.6507E-08 

V 

a 

1.8689E-09 

1.869AE-09 

1.B630E-09 

1.8611E-05 


CALIBRATION OUTPUT PARAMETERS 
TRANS • B*9319E-0l 

RAOHOTC • 3.0296E-09 

0^ RAOCOLC • 1.3627E-05 

H RAOHOTV • 3.5813E-05 

RAOCOLV ■ 1.7151E-05 

OVHOT • 6.1833E-07 

DVCOlt) • -5.8639E-07 
VHOT ■ 2.9987E-05 

VCOLO • 1.3920E-05 

BALANCE OUTPUT PARAMETERS 
RAOHOT ■ 3.7086E-05 

PAOCOLO • 1.5933E-05 

RADHOTG ■ 3.0736E-05 

RADCOLG > 1.3206E-09 

DVHOTB • 2.8268E-09 

OVCOLDB • 2.8266E-09 

VH0T8 • 3.073AE-05 

VCOLOB • 1.3209E-09 


5.3190E-08 6.1731E-08 

1.9^^3E-05 1.9A33E-05 


2.5027E-05 2.A987E-0'5 
2.0867E-05 2.0850E-05 
2.7792E-05 2.777AF-05 
1.2760E-07 1.A319E-07 
2.O0O3E-O5 2.0778E-09 


2.239^E-05 2.2372E-05 
1.8531E-05 1.B519E-05 
2.5556E-05 2.5343E-09 
7.3629E-08 7.8779E-08 
1.859^E-05 1.8379E-05 


6.933<»E-O0 7.6165E-08 

1.9A2<tE-05 1.9A15E-05 


2.<»951E-05 2.^9l7E-05 
2.0833E-05 2.0817E-05 
2.7757E-05 2.77^1E-05 
1.56A2E-07 1.67B7E-07 
2.0755E-05 2.0733E-05 


2.235^E-05 2.2336E-05 
1.8607E-05 1.8596E-05 
2.5531E-05 2.5520E-05 
B.269^E-08 8.578^E-08 
1.8566E-05 1.8553E-05 


e.23^8E-08 9.3129E-08 

1.9A07E-05 1.9392E-05 


2.^805E-O5 2.^8306-05 
2.0801E-05 2.0771E-05 
2.7725E-05 2.7695E-05 
1.7791E-07 1.9^7^E-07 
2.0712E-05 2.067^E-05 


2.2320E-05 2.2291E-05 
1.8585E-05 1.856^6-05 
2.5509E-05 2.5^89E-05 
8.8285E-08 9.2072E-06 
1.B5A1E-05 1.B516E-05 


UYEil 3 

THE SPECTRAL LINE PROFILE IS LORENTZ 


ATMOSPHERIC OUTPUT PARAMETERS 

THE CONCENTRATIONS OF THE INVESTIGATED POLLUTANT CONSIDERED ARE 



D« 

l.OOOOE-07 

2.0000E-07 

3.0000E-07 

A.OOOOE-07 

5.0000E-07 

6.0000E-07 

7.0000E-07 

B.OOOOE-07 

l*0000E-06 

TRANS 

■ 0.3636E-O1 

8.256AE-01 

8.16A3E-01 

8.0838E-01 

8.0123E-01 

7.9A81E-01 

7.8898E-01 

7.8362E-01 

7.7865E-01 

7.6966E-01 

TOTATM 

• 7.90eiE-06 

8.3869E-06 

8.7973E-06 

9.1950E-06 

9.A715E-06 

9.7953E-06 

I.0013E-05 

1.02A9E-05 

1.0A68E*09 

i.0663E-09 

SURFACE 

TOTSUN 

EMISSIVITY - 
• A.6891E-07 

• 980 

A.A607E-07 

A.3250E-07 

^.22A2E-07 

A.1A11E-07 

A.0692E-07 

A.00A9E-07 

3.9^63E“07 

3.8922E-07 

3.79A5E-07 

SURFACE 

TOTSUN 

EMISSIVITY • 

• 2.8135E-06 

• 860 

2.676AE-06 

2.9990E-06 

2.93A5E-06 

2.A8A7E-06 

2.AA19E-06 

2.A029E-06 

2.367BE-06 

2.3353E-06 

2.2767E-06 

SURFACE 

SURFACE 

TOTSURF 

EMISSIVITY - 
TEMPERATURE • 
• ^.90A7E-05 

• 980 

300*000 

A.8A29E-09 

A.7890E-09 

A.7A22E-09 

^•7006E-05 

A.6633E-09 

A.6293E-05 

^.5982E-05 

^.9693E-09 

A.9169E-09 

TOTRAD 

- 5.7A2AE-05 

9.7298E-09 

9.7119E-09 

5.6999E-05 

9.6892E~09 

5.6795E-05 

5.6707E-05 

5.6625E-09 

5.694f9E-05 

5.6A12E-05 

SURFACE 

SURFACE 

TOTSURF 

EMISSIVITY - 
TEMPERATURE • 
• A*2688E*05 

• 980 

296.000 

A*21A7E-09 

4.1681E-05 

^.127AE-05 

A.0912E-05 

A.O5e0E-O5 

A.0292E-05 

A.0021E-r05 

3.9769E-05 

3*931AE-09 

TOTRAD 

• 9»106SE-09 

9«0980E-09 

9.09UE-09 

5.O091E-O5 

S.0798E-05 

5.07S0E-09 

9.0706E-05 

5.066AE-09 

5.0626E-09 

9*0596E-05 

SURFACE 

SURFACE 

TOTSURF 

EMISSIVITY » 
TEMPERATURE • 
• A.AOAZE-09 

• 880 

300.000 

A.3A8AE-09 

A.3003E-09 

A. 2583^-09 

A.2210E-05 

<i.l075F-O5 

A.1570E-05 

<i.l290E-05 

A.1030E-05 

^•0560E-09 

TOTRAD 

• 5.A76AE-09 

9.A.9A7E-09 

9.4395E-05 

5.A272E-09 

5.A166E-05 

5.A071E-05 

5.3985E-05 

5.3906E-09 

9.3833E-05 

5.3700E-05 

SURFACE 

SURFACE 

TOTSURF 

EMISSIVITY « 
TEMPERATURE • 
■ 3.8332E-09 

• 880 

296.000 

3.70A6E-O9 

3.7428E-05 

3.7062E-05 

3.6739E-05 

3.6AA6E-09 

3.6181E-09 

3.5937E-09 

3.5711E-05 

3.5302E-05 

TOTRAD 

• A.905AE-05 

A.09O9E-O5 

A.e820E*05 

<r.8752E-05 

A.869^E-05 

A.86A3E-05 

A.B597E-05 

A.055AE-O5 

A.851AE-05 

A.8^4i2E-05 


INSTRUMENT OUTPUT PARAMETERS 
_ RAOVACB • ^.082AE-05 

^ RAOGASB • 3*3833E-05 

TAUA • 8.2869E-01 

ERRON • 2.6392E-09 


THE CONCENTRATIONS OF THE INVESTIGATED POLLUTANT CONSIDERED ARE 




0. . 

l.OOOOE-07 

2.0000E-07 

3.0000E-07 

4.0000E-07 

5.0000E-07 

6.0000E-07 

7.0000E-07 

8.0000E-07 

l.OOOOE-06 

TRANS 

• 

8.6226E*-01 

8.A226E-01 

0.6Z26E-OI 

8.A226E-01 

8.^Z26E-01 

0.AP26E-OI 

B.^226E-01 

e.^2Z6E-01 

8.A226E-^01 

8.6Z26E-01 

SURFACE 

EMISSIVITY • 

• 980 









SURFACE 

TEMPERATURE • 

300.000 









RAOVACA 

■ 

2.6297E-09 

2.6209E-09 

2.6136E-09 

2.6073E-05 

2.6017E-05 

2.5967E-05 

2.5921E-05 

2.5B79E-05 

2.5839E-05 

2.5768E-05 

RADGASA 

• 

2.178AE-09 

2.1767E-09 

2.1790E-09 

2.1733E-09 

2.I717E-05 

2.1701E-05 

2.1685E-05 

2.1670E-05 

2.l65'iE-09 

2.1625E-05 

RADGASE 

u 

2.07O9E-O9 

2.8691E-09 

2.867AE-09 

Z,8658E-09 

2.86^lE-05 

2.8625E-05 

2.8610E-05 

2.859^E-05 

2.8579E-05 

2.85^9E-05 

DV 

• 

-7.985AE-09 

A.7781E-08 

9.10AAE-08 

1.2663E-07 

1.5662E-07 

1.8231E-07 

2.0A60E-07 

2.2A15E-07 

2.^1^5E-07 

2.7076E-07 

V 

m 

2.1788E-09 

2.17A3E-09 

2.170AE-09 

2.1670E-05 

2.1638E-05 

2.1610F-05 

Z.1583E-05 

2.1558E-05 

2.153<rE-05 

2.1A09E-O5 

SURFACE 

EMISSIVITY - 

• 980 









SURFACE 

TEMPERATURE ■ 

296.000 









RADVACA 

• 

2.3397E*09 

2.3312E-09 

2*3276E-05 

2.32ASE-05 

2.3217E-05 

2.3192E-05 

2.3169E-05 

2.31<t8E-09 

2.3128E-05 

2*3092E-05 

RADGASA 

• 

1.935AE-09 

1.93A9E-09 

1.9336E-09 

U9327E-09 

1.9319E-05 

1.93UE-05 

1.9302E-05 

1.9295E-05 

1.9287E-05 

1.9272E-05 

RADGASE 

a 

2.6278E-09 

2.6269E-09 

2.6260E-05 

2,6292E-05 

2.62<»3E-05 

2.6235F-05 

2.6227E-05 

2.6219E-05 

2.6211E-05 

2.6196E-05 



ov 

■ 

-2«H3^E-09 

2.6095E-08 

4.7206E-08 

6.*t<i^0E-08 

V 

• 

1.9359E-05 

1.9331E-09 

1.9312E-05 

1.9295E-05 

surface 

EMISSIVITY • 

• 880 



SURFACE 

TEMPERATURE ■ 

300*000 



RAOVACA 

■ 

2.5077E-05 

2.A962E-09 

2.4883E-05 

2.4819E-05 

RADGASA 

■ 

2.0772E-09 

2 .07478-09 

2.0723E-05 

2.0701E-05 

RAOGASE 

■ 

2.7697E-05 

2.7672E-09 

2.7648E-09 

2.7625E-05 

OV 

• 

•8.7233E-09 

6.1094E-08 

1.0302E-07 

1.3348E-07 

V 

■ 

2.0777E-05 

2.0717E-09 

2.0672E-05 

2.0634E-05 

SURFACE 

EMISSIVITY • 

.880 



SURFACE 

TEMPERATURE • 

296.000 



RAOVACA 

■ 

2.2^37E-05 

2.2361E-05 

2.2315E-05 

2.2279E-05 

RAOGASA 

■ 

1.8590E-05 

1.8572E-05 

1.8556E-05 

1.85A0E-05 

RAOGASE 

• 

2.551AE-09 

2 .54968-05 

2.9480E-05 

2.5465E-05 

OV 

■ 

-3.8366E-09 

4.1621E-08 

6.3659E-08 

7.7635E-08 

V 

■ 

1.8592E-05 

1.8991E-09 

1.8524E-05 

1.8501E-05 


CAIIBRATION OUTPUT PARAMETERS 
TRANS - 8i9319E-01 

RAOHOTC « 3.0296E-05 

RAOCOLC • l*3627E-05 

RADHOTV - 3.5013E-O5 

RAOCOLV • 1.7151E-05 

OVHOT - 6.1033E-O7 

DVCOlO • -9.8639E-07 
W VHOT « 2.9987E-05 

VCOLO - 1.3920E-05 

balance output PARAMETERS 
RADHOT • 3.7086E-05 

RAOCOLO • l»5933E-09 

RAOHOTG • 3.0736E-09 

RAOCOLG • 1.3206E-05 

OVHOTB • 2.8268E-09 

OVCOlOB • 2.8268E-09 

VHOT8 • 3.073^E-05 

VCOLOB • l,3205E-09 


7.9026E-08 9.1621E-08 

1.9279E-05 1.9265E-05 


2.A76<iE-05 2.^71AE-05 
2.0679E-05 2.0658E-05 
2.7603E-05 2.7582E-05 
1.5750E-07 l.7725E«07 
2.0600E-05 2.0569E-05 


2,22^9E-05 2.2223E-05 
1.B526E-05 1.8512E-09 
2.5^90E-05 2«5^36E-05 
8.7821E-08 9.5818E-08 
1.8^82E-05 1.8A6AE-05 


l,0265E-07 1.12^IE-07 

1.9251E-05 1.9238E-05 


2,A670E-05 2.^629E-05 
2.0638E-05 2.0618E-05 
2.7562E-05 2.75^3E-05 
l,9392E-07 2.082<»E«07 
2.05^1E-05 2.051AE-05 


2.2199E-05 2.2177E-05 
1.8A98E-05 1.8486E-05 
2.5^23E-05 2.5^10E-05 
1.0238E-07 1.0791E-07 
1.8AA7E-05 1.8A32E-05 


1.2H3E-07 1.3607E-07 

1.9226E-05 1.920^E-05 


2.^59lE-05 2.<»522£-05 
2.0599E-05 2.0563E-09 
2.752AE-05 2.7A87E-05 
2.2071E-07 2.^1^lE-07 
2.0<»89E-05 2.0^A2E-05 


2,2156E-05 2.2118E-05 
l.B<»73E-05 1.8^50E-05 
2.5398E-05 2.537^E-05 
1.1267E-07 Ji.20<»8E-07 
1.8^17E-09 1.8390E-05 



LAYER h 

THE SPECTRAL LINE PROFILE IS LORENTZ 


ATMOSPHERIC OUTPUT PARAMETERS 



0« 

l.OOOOE-07 

2.0000E-07 

3.0000E-07 

A.OOOOE-07 

5.0000E-07 

6.0000E-07 

7.0000E-07 

e.OOOOE-07 

l.OOOOE-06 

TRANS 

• 0.15Z4E-O1 

8.0183E-01 

7.9079E-01 

7.81^6E-0l 

7.7336E-01 

7.6620E-01 

7.5976E-01 

7.5300E-O1 

7.^8^6E-01 

7.3868E-01 

TOTATM 

■ 8,67^5E-06 

9.2558E-06 

9.73l^E-06 

1.0132E-05 

1.0^78E-05 

1.0782E-05 

1.1056E-05 

1.1305E-05 

X.153AE-05 

1.1966E-05 

SURFACE 

TOTSUN 

EMISSIVITY • 

■ <t.6215E-07 

• 980 

A.3057E-O7 

A.2A0OE-O7 

^.1A50E-O7 

6.0616E-0.7 

3.9887E-07 

3.9235E-07 

3.86A1E-07 

3.8093E-07 

3.7101E-07 

SURFACE 

TOTSUN 

EMISSIVITY • 

• 2.77Z9E-06 

• 800 

Z.631AE-06 

2.5A68E-06 

2.A875E-06 

Z.A370E-06 

2.393ZE-06 

2.35A1E-06 

2.318AE-06 

Z.Z056E-O6 

2.2261E-06 

SURFACE 

surface 

TOTSURF 

EMISSIVITY • 
TEMPERATURE - 
■ 4.7735E-05 

.980 

300.000 

A.6958E-05 

A.6310E-O5 

^.5777E-05 

6.5307E-05 

A.A891E-05 

9.<»517E-05 

A.A175E-05 

4.3660E-05 

6.3292E-05 

TOTRAD 

• 5.6072E-O5 

5.6653E-05 

5.6A75E-05 

5.6323E-05 

5.6190E-05 

5.6072E-05 

5.5965E-05 

5.5066E-O5 

5.5775E-05 

5.5610E-05 

surface 

SURFACE 

TOTSURF 

EMISSIVITY - 
TEMPERATURE ■ 
■ A.15^5E-05 

• 980 

296.000 

A.0869E-05 

^.0312E-05 

3.98^1E-05 

3.9^32E-05 

3.9070E-05 

3.87A5E-05 

3.0'»A0E-O5 

3.0176E-O5 

3.7679E-05 

TOTRAD 

• 5.0682E-05 

5.0563E-05 

5.0A69E-05 

5.0387E-05 

5.0316E-05 

5.0291E-05 

5.0193E-05 

5.0139E->05 

5.OO80E-O5 

4.9997E-05 

SURFACE 

SURFACE 

TOTSURF 

EMISSIVITY • 
TEMPERATURE ■ 
■ ^.206^E-O5 

.080 

300.000 

A.2167E-05 

A.1592E-05 

^.1106E-05 

A.068^E-05 

^.0310E-05 

3.997^E-05 

3.9668E-05 

3.9305E-O5 

3.087^E-O5 

TOTRAD 

• 5.A311E-05 

5.^056E-05 

5.3872E-05 

5.3725E-05 

5.3590E-O5 

5.3A86E-05 

5.3389E-05 

5.3291E-05 

5.3209E-05 

6.30^7E-05 

SURFACE 

SURFACE 

TOTSURF 

EMISSIVITY • 
TEMPERATURE » 
• 3.7306E-05 

.000 

296.000 

3.6699E-05 

3.6199E-05 

3.5776E-05 

3.5A08E-05 

3.508^E-05 

3.^791E-05 

3.^52^E-05 

3.6278E-05 

3.3839E-09 

TOTRAD 

• A.0753E-O5 

6.8586E-05 

f .0^79E-O5 

A.8395E-05 

A.8323E-05 

9.B259E-09 

6.0ZO1E-O5 

^.81^7E-05 

A.8098E-05 

4.8007E-05 


G\ 


INSTRUMENT OUTPUT PARAMETERS 
RAOVACB - ^.O02^E-O5 

RAD6AS0 - S.3B33E-05 

TAUA • 0.2869E-O1 

ERRON • Z.ea'JZE-O*? 

THE CONCENTRATIONS OF THE INVESTIGATED POLLUTANT CONSIDERED ARE 




0. 

l.OOOOE-07 

2.0000E-07 

3.0000E-07 

^i.OOOOE-07 

5.0000E-07 

6.0000E-07 

7.0000E-07 

8.0000E-07 

l.OOOOE-06 

TRANS 

m 

8.<»226E-01 

8.^226E-01 

8.A226E-0I 

0.6226E-O1 

0.^226E-O1 

8.A226E-01 

B.A226E-01 

8.<»226E-01 

0.9226E-O1 

6.9226E*0l 

SURFACE 

EMISSIVITY • 

• 900 









SURFACE 

TEMPERATURE • 

300,000 









RADVACA 

a 

2.6062E-05 

2.9966E-0$ 

2,5852E-05 

2.5773E-05 

2.570^E-05 

2.566ZE-05 

2.55B6E-05 

2.5535E-05 

Z.5680E-O5 

2.5^03E-05 

RADGASA 

a 

2.15S8E-05 

2.1968E-05 

2.1992E-05 

2.1520E-05 

Z.1^99E-05 

2.1A70E-O5 

2.1A57E-05 

Z.l<i37E-05 

2.1910E-O5 

2.1379E-05 

RADGASE 

a 

2.0512E-O5 

2.8689E-09 

2.8667E-05 

Z.0^A5E-O9 

2.8923E-05 

2.8A02E-05 

2.8382E-05 

2.8362E-05 

2.8392E-05 

2.8306E-05 

DV 

a 

-9.9065E-09 

6.3266E-08 

1.105OE-O7 

1.62A3E-07 

1.98^1E-07 

2.285ZF-07 

2.5415E-07 

2.7628E-07 

2.9563E-07 

3.2792E-07 

V 

a 

2.1593E-05 

2.1533P-05 

2.1S83E-05 

2.1939E-05 

2.U00E-05 

2.136AE-05 

2.1330E-05 

2.1299E-05 

2.1270E-05 

2.1Z15E-05 

SURFACE 

EMISSIVITY • 

• 980 









surface 

TEMPERATURE ■ 

296.000 









RADVACA 

a 

2.3193E-05 

2.3131E-05 

2.3081E-05 

2.3O30E-O5 

2.3001E-05 

2.2967E-09 

2.2937E-05 

2.2909E-05 

2.2803E-O5 

2.2835E-05 

RADGASA 

a 

1.9216E-05 

1.920AE-05 

1.9192E*09 

1.9180E-05 

1.9I60E-O5 

1.9197E-05 

1.9146E-05 

1.9135E-05 

1.9125E-05 

1.9109E-05 

RADGASE 

a 

Z.61A1E-05 

2.6128E-09 

2.6116E-09 

2.610AE-05 

2.6093E-05 

2.6O01E-O5 

2.6070E-05 

2.6060E-05 

2.6069E-05 

2.6029E-05 



DV 

• 

-3.77A0E-09 

3.552AE-0B 

6.^73^E-08 

8.0169E-O8 

1.076AE-07 

1.2A16E-07 

1.3893E-07 

1.5090E-07 

1.6X92E-07 

1.8059E-07 

V 

• 

1.9216E-05 

1.9186E-05 

1.9159E-05 

1.9136E*09 

la911AE-05 

1.9095E-05 

1.9077E-05 

1.9060E-05 

1.9066E-05 

1.901<»E-05 

SURFACE 

EMISSIVITY ■ 

«880 







SURFACE 

TEMPERATURE ■ 

300.000 









radvaca 

■ 

2.A803E-O5 

2.47A7E-05 

2.^652E-05 

2.6575E-05 

2.A509E-05 

2.AA51E-05 

2.^3986-05 

2.A350E-05 

2.A305E-05 

2.422^E-05 

RADCASA 

■ 

2.0610E-05 

2.0581E-05 

2.0552E-05 

2.0526E-05 

2.0500E-05 

2.0A75E-05 

2.0A51E-05 

2.O^20E-O5 

2.0%05E-05 

2.0362E-05 

RAOGASE 

■ 

2.7535E-05 

2.7505E-05 

2.7^77E-05 

2.7650E-05 

2.7^2AE-05 

2.7AO0F-05 

2.7376E-05 

2.7352E-05 

2.7330E-05 

2.7287E-05 

DV 

a 

-1.0501E-08 

7.2720E-08 

1.23A2E-07 

1.6032E-07 

1.8926E-07 

2.1288E-07 

2.3265E-07 

2.4951E-07 

2.6A10E-07 

2.881^6-07 

V 

a 

2.0615E-05 

2.05AAE-05 

2.0A91E-05 

2.0^A5E-05 

2.0^05E-05 

2.0169F-05 

2.0335E-05 

2.0303E-05 

2.0273E-05 

2.0218E-05 

SURFACE 

EMISSIVITY - 

• 880 








SURFACE 

TEMPERATURE • 

296.000 









RAOVACA 

a 

2.2307E-05 

2.2219E-05 

2.2163E-05 

2.2U9E-05 

2.2082E-05 

2.20A9E-09 

2.2019E-05 

2.1991E-05 

2.1965E-05 

2.191BE-05 

RADGASA 

a 

1.8A81E-05 

1.0A6OE-O5 

1.8AA2E-05 

1.8^2^E-09 

1.8A07E-05 

1.8391E-05 

1.8376E-05 

1.8361E-05 

1.83^6E-05 

1.8319E-05 

RAD6ASE 

a 

2.5A05E-05 

2.9385E-05 

2.5366E-05 

2.53A8E-05 

2.5331E-05 

2.5315E-05 

2.5300E-05 

2.5285E-05 

2.5271E-05 

2.5246E-05 

OV 

a 

-A.99A9E-09 

A.7011E-O8 

7.51«6E-08 

9.3638E-08 

1.0775E-07 

1.1917E-07 

1.287AE-07 

1.3692E-07 

1.6<»0^E-07 

1.558AE-07 

V 

a 

1.8A83E-03 

1.0A37E-O5 

I.8A0^E-05 

1.0377E-O5 

1.0353E-O5 

l.«331E-05 

1.8311E-05 

1.0292E-O5 

1.8276E-05 

1.B2A1E-05 

CAtIBRATION OUTPUT PARAMETERS 









TRANS 

a 

8.9319E-01 










RAOHOTC 

a 

3.0296E-05 










RADCOLC 

a 

1.3627E-09 










RAOHOTV 

a 

3.5013E-O5 










RADCOLV 

a 

1.7151E-05 










OVHOT 

a 

6.1833E-07 










OVCOLO 

a 

-9.8639E-07 










VHOT 

a 

2.9907E-O5 










VC 010 

a 

1.3920E-05 










BALANCE 

OUTPUT parameters 









RAOHOT 


3.7O06E-O5 










RADCOLD 


1.5933E-05 










RAOHOTG 


3.0736E-05 










RADCOLG 


1.3206E-05 










OVHOTB 


2.0268E-O9 










OVCOlOB 


2.826BE-09 










VHOTB 


3.073AE-05 










VCOLOB 


1.3209E-05 












APPENDIX C 


Unit Conversion 

The authors wish to note that in the past various units 
have been used to define the line intensities; therefore, the 
following conversion factors may be helpful. 

At standard temperature and pressure condition, STP, 

1 (cm-atm) g,pp = 2.69 x 10^® molecules/cm^ ; 

however, at some temperature T' , 

1 (cm-atm)^, = ^ x 2.69 x molecules/cm^ 
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